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ABSTRACT 


The  sensitivity  of  the  XR-3  roll  behavior  in  calm 
sea  on  sidewall  parameters  and  force  and  moment 
calculations  is  investigated  with  the  Loads  and 
Motions  computer  program  and  compared  with 
experimentally  measured  data.  Propulsion  and  rudder 
subroutines  and  added  mass  computation  are  reviewed 
and  modified.  Recommendations  for  improved  simulation 
of  the  XR-3  roll  behavior  are  given. 
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I.   INTRODUCTION 


A.  Background 

In  the  past  several  years  the  Surface  Effect  Ship  (SES) 
has  received  increased  attention  in  the  United  States  Navy 
and  detailed  investigations  of  its  sea  performance  have  been 
carried  out.  The  two  major  categories  of  SES  are  the  air 
Cushion  Vehicle  (ACV)  and  the  Captured  Air  Bubble  Ship 
(CAB) .  The  CAB-type  of  craft  has  become  of  primary  interest 
to  the  United  States  Navy  and  two  of  the  constructed 
testcrafts  are  the  Bell  Aerospace  Systems  100-B  and  the  XR-3 
of  approximately  100  and  3  tons  of  displacement, 
respectively.  In  1976  the  U.S.  Navy  Surface  Effect  Ship  SES 
100-B  of  the  CAB-type  established  a  world  speed  record  for 
surface  type  ships  of  89.4  knots.  In  April  1976  the  SES 
100-B  launched  a  missile  at  60  knots.  The  launch  was 
successful  and  the  missile  hit  its  target  five  miles  away. 

In  order  to  investigate  the  characteristics  of  the 
CAB-type  SES  under  any  design  and  operating  conditions 
without  the  costly  need  of  modifying  the  actual  craft,  a 
Loads  and  Motions  digital  computer  program  was  developed  by 
Oceanics,  Inc.  [Ref.  1].  In  1970  the  XR-3  testcraft  was 
delivered  to  the  Naval  Postgraduate  School  and  Leo  and 
Boncal  converted  the  1 00-B  Loads  and  Motions  program  to 
represent  the  XR-3.  Since  there  were  substantial  design 
differences  between  these  two  ships  programming 
modifications  were  required  in  certain  subroutines. 


B.  Objectives 

The  purpose  of  this  thesis  is  to  investigate  one  aspect 
of  safe  maneuverability  of  these  high  speed  ships,  i.e. 
their  roTl  behavior  in  a  turn  maneuver.  To  fulfill  this 
objective  the  sensitivity  of  the  simulated  XR-3  roll 
behavior  in  turn  maneuvers  as  it  is  represented  by  the  Loads 
and  Motions  Program  is  investigated  for  its  dependence  on 
changes  in  sidewall  characteristics.  The  effect  of 
modifications  in  the  added  mass  computation,  sidewall 
parameters  and  force  calculations  as  well  as  propulsion  and 
rudder  input  parameters  are  investigated.  An  evaluation  of 
the  computed  performance  is  obtained  by  comparison  with 
experimentally  measured  steady  state  roll  conditions. 
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II.   INITIAL  REMARKS 


A.   PROBLEM  OF  REPEATABILITY 


It  has  been  this  author's  experience  as  well  as  all 
previous  investigators  at  the  Naval  Postgraduate  School 
using  the  XR-3  Loads  and  Motions  Program  that  there  have 
been  no  observed  problems  related  to  repeating  the  results 
for  a  given  input  condition.  However,  during  the  author's 
first  studies  of  the  Loads  and  Motions  Program  as  listed  in 
Menzel's  thesis  [Ref.  2]  an  attempt  was  made  to  duplicate 
some  of  the  simulation  runs  given  in  that  reference  work.  It 
was  found  necessary  as  a  first  step  to  restore  the 
simulation  program  on  the  IBM  360/67  computer  from  a  card 
source  onto  a  disk. l  Using  the  input  conditions  stated  in 
Ref.  2  it  was  impossible  to  obtain  identical  time 
histories,  e.g.  for  turns.  Each  time  subroutine  SIDEWALL  was 
read  into  the  computer  together  with  the  data  input  deck, 
thereby  overriding  the  stored  program,  an  error  message  in 
the  printed  output  indicated  that  a  division  by  zero 
occurred  in  this  subroutine.  The  error  message  was  not 
generated  when  the  stored  program  which  contained  an 
identical  version  of  subroutine  SIDEWALL  was  called.  A 
closer  look  revealed  that  in  the  SIDEWALL- version  given  in 
Ref.  2   the   variable   PBAR   (plenum  pressure)   was   neither 


1  It  should  also  be  noted  that  results  obtained  in   Ref.  2 

used  the  IBM  360/67  Release  20.6  system  while  the  results 

obtained  in  this  work  used  the  21. 8B   Release  with   HASP 
installed. 
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defined  nor  transferred  by  a  COMMON-statement  and  therefore 
a  default  value  of  zero  was  used  during  the  computation  of 
PBHEAD  (see  Appendix  C,  SDWL  0440)  .  The  missing  statement 

PBAR  =  PB  -  PINF 

was  inserted  to  subroutine  SIDEWALL.  With  this  correction 
the  new  results  were  slightly  closer  to  those  given  in 
Ref.  2. 

The  table  given  below  compares  output  values  for  roll 
and  pitch  angle  in  a  turn  at  20  kn  with  constant  thrust  and 
a  35  degrees  rudder  step  input  under  calm  sea  conditions. 

TABLE_I 
20  kn  turn  at  35  deg  rudder  angle 

without    with  relative 

PBAR-card  PBAR-card   Ref.  2   deviation 

Roll  angle  (deg) 


first  peak 

7.44 

7.39 

7.2 

2.6 

% 

avrg  at  20  sec 

.  2.65 

2.77 

3.  1 

11.9 

7, 

itch  angle  (deg) 

first  peak 

1  .34 

1  .33 

1.2 

9.3 

0 

avrg  at  20  sec 

0.99 

0.94 

0.8 

14.9 

% 

Although  the  differences  between  columns  2  and  3  in  TABLE  I 
might  be  considered  to  be  small  as  far  as  magnitudes  are 
concerned  ,  the  relative  deviation  is  up  to  15  %  as  shown  in 
column  4.  But  what  is  more  important  is  the  fact  that  the 
roll  and  pitch  angle  responses  for  the  uncorrected  (plots  1 
and  2  in  Appendix  A)  as  well  as  for  the  corrected  program 
(plots  3  and*  4)  show  an  unstable  craft  behavior  for  t>40  sec 
where  the  pitch  angle  increases  rapidly  while  the  roll  angle 
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decreases  rapidly  approaching  the  zero  degree  value.  This 
unstable  condition  is  probably  due  to  the  35  degrees  rudder 
step  input  at  20  kn,  as  it  is  used  in  Menzel's  study,  which 
is  physically  unreal  since  it  could  never  be  introduced  to 
the  XR-3  testcraft  in  an  actual  run  at  that  speed. 
Therefore  rudder  angles  of  up  to  15  degrees  introduced  at  a 
rate  of  5  deg/sec  will  be  used  throughout  this  work,  as  it 
is  also  done  on  the  testcraft. 

Also,  regarding  plots  1-4  again  and  not  considering 
whether  they  really  reflect  the  actual  craft  behavior,  it 
should  be  noted  that  a  short  time  history  (  t<20  sec  in 
Ref.  2)  could  possibly  result  in  wrong  conclusions  because 
the  unstable  condition  is  not  evident  at  that  instant  of 
time.  Therefore,  time  histories  of  up  to  40  or  50  seconds 
will  be  shown  throughout  this  study. 
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B.   STEADY  STATE  CONDITIONS 


Due  to  the  change  in  subroutine  SIDEWALL  as  mentioned  in 
the  preceding  section  the  steady  state  conditions  for 
straight  runs  in  calm  water  have  been  reestablished  for 
various  speeds  and  are  listed  in  Table  II. 

Table  II 
Steadv  State  Conditions 


Speed 

Draft 

Pitch 
angle 

Plenum 
pressure 

Thrust 

(kn) 

(in) 

(deg) 

(psf) 

(lb) 

10.0 

8.  17 

1.62 

23.93 

400.61 

12.5 

7.03 

1.11 

24.84 

349.42 

15.0 

6.74 

0.84 

2U  .84 

335.71 

17.5 

6.41 

0.63 

24.84 

346  .45 

20.0 

6.12 

0.48 

24.84 

373.44 

22.5 

5.87 

0.36 

24.  84 

411 .53 

25.0 

5.66 

0.29 

24  .84 

458.62 

27.5 

5.48 

0.25 

24.84 

513.31 

30.0 

5.34 

0.26 

24.84 

574  .43 

These  steady  state  conditions  have  been  established  by 
first  executing  the  XR-3  Loads  and  Motions  Program  with 
constant  speed  for  40  seconds  and  then  repeating  the  run 
keeping  the  evaluated  thrust  constant.  A  comparison  of 
these  steady  state  values  with  those  given  in  Ref.  2  shows 
differences  especially  in  the  lower  and  upper  speed  range  of 
again  up  to  15  %  the  cause  of  which  could  not  be  suspected 
anywhere  else  since  the  simulation  program  has  not  been 
changed  since  last  used  by  Menzel. 
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C.       SIGN    CONVENTIONS 


The  sign  conventions  used  in  the  XR-3  Loads  and  Motions 
Program  are  identical  to  those  used  in  Report  71-84  by 
Oceanics  Incorporated  [Ref.  1].  A  right  handed  coordinate 
system  is  applied  to  the  craft  with  the  positive  X-axis 
being  measured  forward  and  the  lateral  Y-axis  being  measured 
positive  to  starboard.  The  vertical  Z-axis  is  measured 
positive  downward.  Identical  signs  are  used  for  respective 
forces   along    those   axes. 

The  signs  of  the  angles  are  defined  in  the  following 
manner    : 

pitch  angle    being   positive   upward 

(boat    noses   up) 

roll  angle  being  positive  to  starboard 

(boat  rolls  to  starboard) 

yaw   angle  being  positive  to  starboard 

(boat  turns  to  starboard) 

rudder  angle  being  positive  to  port 

(right  rudder,  boat  turns  to  port)  . 

Zero  pitch  and  roll  angle  are  referenced  to  the  X-Y 
plane  parallel  to  the  water  surface. 
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III.       INVESTIGATION    OF    SIDEWALL    EFFECTS 


A.       ADDED    MASS    EFFECT 


1 .       Theory 

A  basic  element  in  representing  the  hydrodynamic 
forces  by  use  of  slender  body  theory  are  the  two-dimensional 
sectional  values  of  lateral  added  mass  A22  of  the  sidewalls, 
which  are  necessary  for  the  determination  of  the  lateral 
forces,  as  well  as  the  two-dimensional  sectional  vertical 
added  mass  A33  which  is  also  used  in  determining  the  roll 
moment.  These  two-dimensional  added  mass  values  are  given  in 
the  'Surface  Effect  Ships  Aero/Hydrodynamics  Technology 
Design  Manual'  [Ref.  3  ]  as 


and 


A22  =  C  *RHO*PI*D*D/2.0 
h 


A33  =  C  *RHO*B*3/S.O 
v 


where  C   is  the  lateral  added  mass  coefficient 
h 

C   is  the  vertical  added  mass  coefficient 
v 

RHO  is  the  specific  density  of  sea  water 

D  is  the  local  draft 

and    B  is  the  width  at  the  local  sidewall  waterline. 
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The  C   and  C   values  originally   selected   in   Ref.  1    were 
h      v 

those   corresponding   to   high  speed  with  C  =0.4  and  C  =1.0. 

h  v 

But  in  order  to  account  for  variations   in   sidewall   shape, 

and  their  influence  on  the  effective  lateral  added  mass,  the 

results  of  the  research  work  described  in  Ref.  3  led   to   an 

average   value   of  C  =0.8.  This  new  coefficient  improved  the 

h 

agreement   between   theoretical   and   experimental   data   as 

stated   in   Ref.  3  .   In   the  present  XR-3  Loads  and  Motions 
Program  a  value  of  C  =0.4  is  used  in  accordance  with   Ref.  1 

h 

and   the   effect   of   letting  C  =0.8  will  be  investigated  in 

h 

later  parts  of  this  thesis. 

Considering  the  equation  for  the  vertical  added  mass 
A33,  a  major  difference  has  been  found  between  that  one 
given  in  Ref.  3  and  the  one  used  in  the  XR-3  Loads  and 
Motions  Program  (e.g.  Ref.  2),  where  in  accordance  with 
Ref.  1  (1971)  there  was  an  additional  Pi-factor  in  the  A33s 
equation.  The  letter  ' s'  indicates  that  the  computation  of 
the  added  mass  is  done  at  the  craft's  stern.  Since  Ref.  3 
(1976)  reflects  the  experimental  and  theoretical  work  being 
done  based  on  Ref.  1 ,  the  A33s  equations  in  subroutines 
SIDEWALL  and  SIDETAB  have  now  been  changed  to  the  new  form 
given  above. 
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2 •   New  Steady  State  Conditions 

The  change  in  the  simulated  XR-3  performance  due  to 
the  reformulated  A33s  computation  can  be  observed  from 
TABLE  III  giving  the  new  steady  state  conditions  for 
straight  runs  in  calm  sea  for  various  speeds  which  are 
identical  for  both  (0.4  and  0.8)  lateral  added  mass 
coefficients. 

TABLE  III 
Steady  state  conditions 
(new  ,A33s^ 

Speed     Draft 

(kn)      (in) 

15.0  6.80 

17.5  6.46 

20.0  6.17 

22.5  5.92 

25.0  5.71 

27.5  5.53 

30.0  5.40 

Comparing  these  steady  state  values  with  those  previously 
given  in  TABLE  II  in  Sect.  II. 3  one  finds  in  general  for  all 
speeds  that 

-  draft  has  increased  by  about  0.05  in  (0.8  %) 

-  pitch  angle  has  increased  by  0.02  deg  (2  to  7  %) 

-  thrust  has  slightly  increased  by  less  than  0.5  %    . 


Pitch 

Plenum 

Thr usr 

angle 

pressure 

(deg) 

(psf) 

(lb) 

0.86 

24.84 

336.30 

0.65 

24.84 

347.23 

0.49 

24.  84 

374. 38 

0.37 

24.  84 

412. 62 

0.30 

24.  84 

459. 90 

0.27 

24.  84 

514. 89 

0.28 

24.84 

576.66 
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3 .   Effect  ia  a^Turn^maneuyer 

The  effect  of  the  change  in  the  formulation  of  the 
vertical  added  mass  in  a  simulated  turn  maneuver  has  been 
studied  next.  After  a  5  sec  straight  run  at  20  kn  in  calm 
sea  a  port  rudder  deflection  was  introduced  at  a  rate  of 
5  deg/sec  and  then  kept  fixed  at  15  degrees  resulting  in  a 
turn  to  port.  The  deadrise  forces  are  computed  at  the 
transom.  In  TABLE  IV  are  shown  the  steady  state  conditions 
and  the  roll  moments  contributed  from  the  bow  seal  (FKBS) , 
stern  seal  (FKSS) ,  sidewalls  (FKSW) ,  rudder  (FKRUD) , 
propulsion  system  (FKP) ,  aerodynamics  (FKAED)  and  plenum 
pressure  (ABPB*PHI*Z) . 

TABLE_IV 
Steady  state  conditions  in_20._kn.turn 
\   5  deg  port  rudder  angle 

deviation 

pitch  angle  (deg)     0.48  0.52  8.3  % 

roll  angle  (deg)     2.0  1.98  -  1.0  % 

draft  (in)      5.86  5.94  1.4  % 

speed  (kn)  19.25  19.27  0.1  % 

moments: 

FKBS  -391. 1  -386.3  -  1.2  % 

FKSS  -   3.6  -3.6  0.0  % 

FKSW  191.0  202.0  5.8  % 

FKRud  -  97.5  -110.0  -12.8  % 

FKP  0.2  0.2  0.0  % 

FKAed  -  54.4  -  53.5  1.7  % 

ABPB*PHI*Z  355.4  351.9  1.0  % 


old    i 

i33s 

0. 

.48 

2 

.0 

5 

.86 

19 

.25 

new 

A3  3  s 

0 

.52 

1 

.98 

5 

.94 

19 

.27 
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From  TABLE  IV  it  can  be  seen  that,  due  to  the 
modified  A33s  computation,  the  counteracting  roll  moments 
from  the  sidewall  (FKSW)  and  rudder  (FKRUD)  changed  by  about 
6  and  13  %,  respectively,  resulting  in  only  a  little  effect 
in  generating  a  tendency  toward  a  smaller  roll  angle  which 
is  directed  out  of  the  turn.  Other  major  contributions  to 
the  roll  angle  are  from  the  bow  seal  (FKBS  directed  into  the 
turn)  and  from  the  plenum  pressure  (A3P3*PHI*Z)  directed  out 
of  the  turn.  The  corresponding  plots  (plots  5-8  for  old 
A33s,  plots  9-12  for  new  A33s)  show  that  with  the  new 
equation  for  A33s  the  pitch  and  roll  angle  response  curves 
are  less  damped  than  they  were  using  the  former  equation  for 
vertical  added  mass  calculation  as  would  be  expected. 

It  should  be  noted  that  the  roll  moment  doe  to 
rudder  is  fairly  large  and  changes  quite  significantly  while 
the  moment  contributed  by  the  propulsion  is  very  small  and 
constant.  The  cause  ol  this  will  be  investigated  in 
chapter  IV. 
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B.   DEADRISE  ANGLE 


1 .   Background 

Surface  Effect  Ships  (SES)  of  the  captured  air 
bubble  type  (CAB)  have  rigid  sidewalls  immersed  into  the 
water,  thereby  -  together  with  the  bow  and  stern  seal  - 
capturing  the  air  in  the  planum  chamber  and  reducing  air 
leakage.  Since  the  sidewalls  of  the  XR-3  are  not  uniform  in 
cross-sectional  shape  throughout  their  length  but  are  curved 
on  the  outboard  side  near  the  bow  similar  to  a  single  hull 
ship,  the  expression  "sidehulls"  would  be  more  appropriate, 
as  can  be  found  in  modern  literature. 

The  importance  of  the  correct  understanding  of  the 
effect  of  the  sidewalls  on  the  craft's  performance  is 
readily  seen  from  a  report  on  SES  Research  with  the  XR-1 
testcraft  [Ref.  5],  As  reported,  the  first  test  series  were 
carried  out  successfully  with  a  45  degrees  sidewall  deadrise 
angle,  which  is  the  angle  between  the  horizontal  plane  and 
the  ship's  outer  sidewall  surface.  In  1964,  after  a 
modification  of  the  craft  to  60  degrees  deadrise  angle  and 
articulating  seals,  an  unstable  roll  condition  occurred 
during  a  turn  at  35-37  knots,  a  maneuver  performed  many 
times  before.  The  testcraft  heeled  out  of  the  turn  ,  nosed 
down  due  to  the  retraction  of  the  bow  seal,  continued  its 
outward  roll  motion  aQd  then  flipped  over.  The  report, 
however,  does  not  mention  the  rudder  action  generating  the 
turn.  After  this  accident,  the  deadrise  angle  had  been 
changed  back  to  45  degrees.  From  this  experience  the 
surface   effect   ship's  roll  stability  should  be  expected  to 


21 


be  sensitive  to  the  deadrise  angle  and  an  investigation  of 
this  sensitivity  is  reported  in  the  following  sections. 

2 •   Forces  due  to  Deadrise  Angle 

The  deadrise  angle,  as  pointed  out  by  Ref.  3  and  4, 
has  a  major  effect  on  the  craft's  roll  motion  and  therefore 
is  an  important  design  consideration.  During  a  turn  the 
craft  is  desired  to  heel  inboard,  thereby  minimizing  the 
apparent  transverse  acceleration  (a  coordinated  turn).  For 
typical  5ES  designs,  it  is  not  possible  to  achieve  perfect 
coordination,  but  it  is  possible  to  achieve  nearly  flat 
turns  if  the  craft  is  designed  properly.  The  principal 
factors  affecting  roll  stability  at  speed  are 

*  sidewall  geometry 

*  seal  characteristics 

*  vertical  location  of  the  center  of  gravity. 

The  forces  acting  at  the  center  of  gravity  station  of  the 
starboard  sidewall  in  a  port  turn  are  shown  in  Figure  1. 
The  centrifugal  force  generated  acts  away  from  the  center  of 
the  turn.  This  force  must  be  counteracted  by  hydrodynamic 
forces  wnich  are  created  by  yawing  the  ship  into  the  turn. 
To  an  observer  aboard  the  craft,  it  would  appear  as  if  the 
craft  is  sideslipping.  Wave  buildup  on  the  outboard  sides  of 
the  sidewalls  increases  the  pressure  there,  while  only  small 
pressure  changes  occur  on  the  inboard  side  of  the  opposite 
sidewall  generating  a  small  (negligible)  vertical  force. 
The  principal  force  arises  from  the  outboard  side  of  the 
sidewall.  The  direction  of  this  force  depends  on  the  slope 
of  the  deadrise  angle.  If  the  resultant  force  passes  above 
the  center  of  gravity  (solid  line  in   Fig.  1)    a   restoring 
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moment  results,  tending  to  heel  the  craft  into  the  turn 
thereby  improving  its  roll  stability.  If  the  deadrise  angle 
is  chosen  to  be  larger  (e.g.  70  deg,  shown  dashed),  the 
resultant  force  passes  below  the  center  of  gravity, 
resulting  in  a  moment  tending  to  heel  the  craft  out  of  the 
turn. 

Thus,  once  the  SES  geometry  and  the  exact  vertical 
location  of  the  center  of  gravity  are  known,  an  approximate 
check  for  stability  in  a  turn  can  be  made. 


Figure  1  -   GENERAL  ROLL  STABILITY  CHECK  IN  A  PORT  TURN 
(C.F.  CENTRIFUGAL  FORCE,  C. G.  CENTER  OF  GRAVITY) 
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3 .   XR-3  Sidewall  Geometry 

From  the  construction  data  entered  into  the  sidewall 
table  (subroutine  SIDSTAB)  Figure  2  has  been  drawn  showing 
the  cross  sections  at  several  selected  stations  of  the 
starboard  sidewall.  Since  in  most  simulation  runs  (which 
will  be  discussed  later)  draft  was  less  than  7  inches,  the 
deadrise  angles  for  all  stations  immersing  into  the  water 
have  been  calculated  by  straight  line  approximation  and  are 
given  in  TABLE  V.  The  stations  are  counted  from  bow 
(station  0)  to  stern  (station  28)  and  are  not  equally  spaced 
for  station  numbers  less  than  11.  The  center  of  gravity  has 
been  located  experimentally  at  station  14.  The  average 
deadrise  angle  of  all  stations  listed  is  64.3  degrees. 

TABLE_7 

Deadrise,  angles £dec[) fQJZ_§ll_§  Nations 

Straight    line    approximation    for    7,f draft 


Station 

Deadrise  an 

gle 

Station 

Deadrise  angle 

5 

67.6 

17 

62.1 

6 

70.8 

18 

64  .0 

7 

57.9 

19 

65.3 

8 

54.  4 

20 

66  .7 

9 

53.7 

21 

58.  1 

10 

53.7 

22 

6.96 

11 

53.7 

23 

70.3 

12 

55.6 

24 

71.1 

13 

57.3 

25 

72.4 

14 

58.5 

26 

73  .8 

15 

59.7 

27 

76.9 

16 

60.9 

28 

76.8 
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Figure  2  -   XR-3  STARBOARD  SIDEWALL  CROSS  SECTIONS 
(STATIONS  AS  INDICATED,  UNITS  IN  FT.) 
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4  .   Eff  e  ct_o  f  _D  eadri  se_  A  ng  le 

From  the  previous  discussion  the  surface  effect 
ship's  stability  behavior  is  expected  to  be  sensitive  to  the 
sidewalls'  deadrise  angle.  In  the  XR-3  Loads  and  Motions 
Program  listed  in  Ref.  2  the  deadrise  angle  at  the  transom 
(78.8  deg)  had  been  used  for  the  computation  of  the  deadrise 
forces.  This  calculation  was  implemented  into  subroutine 
SIDEWALL  using  the  following  statements: 

FYH  (J)=-A22S*U*  (V+XSS*R-ZS*P) 

• 

CTNDR=0.0 

IF  (DSS.IE.0.0)     GO    10    22 
CTNDR=  (ES-3B(  1)  )  /DSS 
IF  (THETA.LT.0.0)     CTNDR=0 . 33 39  1 
22    CONTINUE 

FZHOLD  (J)=FZH  (J) 

FZHDRP  (J)  =PS1*FYH  (J)  *CT  NDR*PR0M01 

FZH  (J)  =FZH(J)  +FZHDRP(J) 


FK=  {FZH  (2) -FZH  (1)  )  *Y 3W  +  FKD-FY*ZS 

The  projection  of  the  lateral  force  FYH (J)  at  each  sidewall 
(J=1  for  Fort,  J  =  2  for  starboard  side)  produced  the 
respective  component  of  vertical  force  FZH,  where  FZHDRP  is 
the  projected  force  and  CTNDR  is  the  cotangent  of  the 
deadriseN angle.  The  roll  angle  had  been  already  taken  into 
account  in  the  computation  of  DSS  (draft  at  the  stern)  .  The 
factor  PRCM01  simply  provided  a  means  to  aroitrarily  change 
the  vertical  projected  force  for  studies  undertaken  in 
Ref.  2.   The  factor  PM1  introduced  rhe  necessary  sign  change 
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in  the  sidewall  force  computation  being  dependent  on  the 
craft's  side  (FM1  =  -1  for  port,  PM1=+1  for  starboard  side). 
The  roll  moment  FK  was  partly  determined  from  the  vertical 
forces  of  both  sidewalls. 

Comparing  the  deadrise  angles  given  in  TABLE  V  it  is 
seen  that  the  transom  deadrise  angle  of  78.8  degrees  is 
unique  and  the  largest  along  each  sidewall  and  no  reason 
could  be  found  why  just  the  transom  deadrise  angle  should  be 
used  for  the  vertical  force  computation.  This  fact  rather 
came  up  when  Leo  and  Boncal  (Ref.  8)  scaled  down  the 
simulation  program  of  the  100-B  testcraft  to  create  the  XR-3 
Loads  and  Motions  Program.  For  the  100-B  testcraft  the 
sidewall  cross  sections  are  uniformly  shaped  throughout  most 
of  its  length.  This  uniformity  does  not  occur  for  the  XR-3, 
as  shown  in  Figure  2.  Therefore  it  has  been  supposed  that 
another  deadrise  angle  which  is  more  representative  for  all 
angles  existing  along  the  XR-3  sidewalls,  e.g.  the  deadrise 
angle  at  the  center  of  gravity  location  (station  14),  could 
possibly  be  more  appropriate  for  the  vertical  force 
computation.  To  investigate  the  effect  of  this  supposition, 
turn  maneuvers  in  calm  sea  at  20  knots  and  various  rudder 
angles  have  been  simulated  for  both  deadrise  angles  and 
lateral  added  mass  coefficients. 

TABLE, VI 

Steady  State  Roil  Angle  (deq)  at  20  kn 
for  various  Rudder  Angles £decQ_ 

Deadrise  force  computation  at      Experimental 
center  of  gravity      testcraft 

data  (Ref. 6) 


0.09 
0.28 
0.58 
1.36 


Rudder 

Deadris 

angle 

transo  m 

C  =0.  4 
h 

5 

0.52 

9 

1.07 

12 

1.50 

15 

1.98 

C  =0.4 

h 

C  =0.8 
h 

0.32 

0.29 

0.78 

0.72 

1.  15 

1.11 

1  .  53 

unstable 
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Table  VI  shows  the  steady  state  values  for  the  roll 
angle  which  are  compared  with  those  measured  experimentally 
in  1974  and  documented  in  Ref.  6.  Since  the  XR-3  Loads  and 
Motions  Program  under  investigation  represents  the  XR-3 
craft  configuration  as  it  existed  in  1974,  Ref.  6  may  serve 
to  check  whether  the  results  produced  by  certain  changes  in 
the  simulation  program  are  favorable  or  not.  It  is  not 
anticipated  to  exactly  match  the  results  obtained  from  the 
simulated  runs  to  the  measured  values  since  for  the 
testcraft  data  the  following  precisions  for  the  measuring 
devices  are  stated  in  Ref.  6  :  pitch  and  roll  angle  ±0.5°, 
rudder  angle  ±1.0°  and  speed  ±1  kn.  From  TABLE  VI  it  is 
obvious   that   the   steady  state  values  for  roll  angle  using 

C  =0.4  as  the  lateral  added  mass  coefficient  and  the   center 
h 

of  gravity  deadrise  angle  for  the  vertical  force  computation 

are  closer  by  25  to  40  %   to  the  measured  angles  than  are  the 

corresponding   roll  angles  considering  the  transom  deadrise 

angle . 

Using   the   center  of  gravity  deadrise  angle  and  the 

lateral   added   mass   coefficient   C  =0.8   as   suggested   in 

h 

Ref.  3,   the   agreement  in  steady  state  roll  angles  could  be 

improved  by  another  3   to   1 0  %      for   rudder   angles   up   to 

12  degrees.   This   improvement   has   also   been   reported  in 

Ref.  3.   Comparing   the   roll   angle   responses   for    turns 

generated   by   a  12  degrees  rudder  angle  to  port  with  C  =0.4 

h 

(plots  21  and  22)  and  C  =0.8  (plots  23    and  24)  it   is   found 

h 

that  the  smaller  lateral  added  mass  coefficient   generates  a 

little  higher  peak  roll  angle  (1.47°)  than  does  the  larger 
coefficient  (1.27°) .  Both  responses  show  the  same  number  of 
oscillations  (eleven  cycles)  until  they  die  out.  The  larger 
lateral   added   mass   coefficient  generates  a  small  negative 
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roll  angle  shortly  after  the  port  rudder  motion  has  been 
introduced.  However,  for  a  15  degrees  port  rudder  angle  the 
simulation  program  showed  an  unstable  response  with 
increasing  amplitudes  of  oscillation  for  pitch  and  roll 
angle  (see  plots  17  and  18).  From  this  the  lateral  added 
mass  coefficient  being  0.8  does  not  seem  to  be  appropriate 
for  the  XR-3  simulation  if  both  sidewalls  are  considered  for 
the  vertical  deadrise  force  computation.  The  effect  of 
using  both  lateral  added  mass  coefficients  will  be 
considered  again  in  Sections  III.D  and  E. 

From  TA3LS  711  it  can  be  seen  that  the  roll   moments 

contributed    by    the   bow   seal   (FKBS),   plenum   pressure 

(ABPB*PHI*Z)  ,  rudder  (FOOD)  and  sidewalls  (FKSW)  change  the 

most.  Comparing  the  effect  due  to  the  change  from  transom  to 

center  of  gravity  deadrise  angle  it  is  found  that 

-  bow  seal  inward  effect  decreases 

-  plenum  pressure  outward  effect  decreases 

-  rudder  inward  effect  decreases 

-  sidewall  outward  effect  increases 

=  outward  roll  angle  decreases. 

The  time  histories  for  simulated  turn  maneuvers  with 

15  degrees   port  rudder   angle   using  center  of  gravity  and 

transom  deadrise  angle  (TA3LE  711)  are  shown  in   Appendix   A 

as   plots   13-16  and   9-12,   respectively.   Comparing  these 

plots  the  effect  of  the  change  from   transom   to   center   of 

gravity  deadrise  angle  and  C  =0.4  can  be  observed  as 

h 


29 


-  pitch  angle  response  being  more  damped 

-  roll  angle  response  being  more  damped 

-  pitch  and  roll  rate  responses  being  more  damped 
and  having  smaller  peak  values. 


TA3LE_VII 
Steady  State ^Values  for  selected  .Deadrise^Angles 


20 

kn  Turnx  calm  5 

eai_15_deg_ 

_port  _  rudder 

_anc[le 

Deadrise 

An 

gle 

Transom 

C.  of  G. 

arbitrar 

ily  chosen 

(de 

grees) 

78.  8 

58.5 

31.6 

28.7 

Pitch  an 

gle 

(deg) 

0.52 

0.45 

0.52 

0  .64** 

Roll  ang 

16 

(deg) 

1.98 

1.53 

0.31* 

-  0.33* 

Draft 

(in) 

5.  94 

5.81 

6.07 

6.52 

Speed 

(kn) 

19.27 

19.37 

19.1 

18.7 

FYSW 

-824.  0 

-821 .9 

-810.0 

-801  .0 

FYRUD 

40.0 

32.7 

39.5 

55.8 

FYAED 

-  19.  5 

-  20.0 

-  19.1 

-  17.6 

R*V*AM 

803.  5 

809.2 

739.6 

762.8 

FKES 

-386.3 

-336.4 

-  68.4 

71  .6 

FKSS 

-   3.6 

-   2.7 

-   0.5 

0.5 

FKSW 

202.0 

212.  1 

175.0 

184.0 

FKRUD 

-110.0 

-  39  .9 

-108.4 

-153.0 

FKAED 

-  53.5 

-  56.0 

-  52.4 

-  46.2 

FKP 

0.  19 

0.14 

0.03 

-   0.03 

ABPB*PHI 

*Z 

351.  9 

273.5 

54.6 

-  57.2 

FYH  (P/S) 

3 9/- 147 

-45/-126 

-78/-95 

-107/-88 

FYD 

-637.  8 

-651.4 

-637.3 

-607.0 

Note  :   *  =  still  decreasing 

**  =  still  increasing,  not  quite  steady  state 
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The  reason  for  the  sidewall  roll  moment  and  thereby 
the  net  roll  effect  not  changing  more  pronounced  (TABLE  VII) 
is  that  not  cnly  the  lateral  force  FYH  must  be  considered  in 
order  to  determine  the  force  F  effecting  roll  stability 
(Sect.  III.B.2,  Fig.  1)  but  also  the  lateral  drag  force  FYD 
has  to  be  taken  into  account.  Considering  the  XR-3  geometry 
using  the  center  of  gravity  deadrise  angle,  as  it  is  shown 
in  solid  lines  in  Fig.  3,  the  lateral  force  FYH  can  be  seen 
to  generate  the  vertical  projected  force  FZHDRP  and  the 
force  F,  which  is  directed  well  above  the  center  of  gravity. 
But  since  FYD  has  to  be  added  vectorially  to  the  force  F  and 
its  strength  being  several  times  larger  than  FYH  (TABLE  VII) 
the  new  resultant  force  F'  determining  the  craft's  roll 
behavior  is  directed  well  below  the  center  of  gravity 
(although  F  being  directed  above  it).  Thereby  an  outward 
roll  angle  is  generated. 


CG 


FYD 

FZH 
A 

F'          ~~"  —   — .  . ,  ^                 ^^ 

FYH 

WL 


Figure  3  -   DEADRISE  AND  DRAG  FORCES  IN  A  PORT  TURN 
(CG=  CENTER  OF  GRAVITY,  WL=  WATER  LINE) 
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Also  from  Fig.  3  the  effect  of  considering  the 
center  of  gravity  versus  the  transom  deadrise  angle  on  the 
force  generating  the  inward  or  outward  roll  moment  can  be 
seen.  The  center  of  gravity  deadrise  angle  and  its 
corresponding  forces  are  shown  in  solid  or  dashed  lines  with 
F  being  directed  above  the  center  of  gravity.  The  transom 
deadrise  angle,  shown  in  dotted  lines,  generates  a  force  F 
being  directed  below  the  center  of  gravity  producing  an 
outward  roll  moment  which  gets  reinforced  by  adding 
vectorially  FYD.  From  this  the  general  stability  check 
yields  a  smaller  outward  roll  moment  and  angle  using  the 
center  of  gravity  versus  the  transom  deadrise  angle.  This  is 
in  agreement  with  the  results  given  in  TABLE  VI. 

Next,  some  simulation  runs  with  arbitrary  deadrise 
angles  chosen  such  that  small  positive  and  negative  steady 
state  roll  angles  result  are  presented  in  order  to,  first, 
find  out  which  deadrise  angle  would  generate  the  moments 
necessary  to  simulate  flat  turns,  and,  second,  verify  the 
trend  of  the  change  in  the  respective  moments  which  has  been 
observed  for  the  change  from  transom  to  center  of  gravity 
deadrise  angle  in  TABLE  VII.  The  value  of  the  lateral  added 
mass  coefficient  in  these  runs  was  0.4  as  it  was  in  all  runs 
listed  in  TABLE  VII.  Considering  the  roll  moments  over  the 
range  of  deadrise  angles  from  78.8  to  28.7  degrees,  the 
moments  due  to  the  bow  seal,  stern  seal,  propulsion  and 
plenum  pressure  change  consistantly  in  magnitude  while  those 
due  to  sidewalls,  rudder  and  aerodynamics  reach  some  extreme 
value  and  then  increase  again.  Regarding  the  steady  state 
roll  angles  listed  in  TA3LE  VII,  it  can  be  seen  that  the 
arbitrary  chosen  deadrise  angles  happen  to  be  nearly 
symmetric  about  the  deadrise  angle  that  would  generate  zero 
roll  angle  (flat  turn)  under  these  simulation  conditions. 
Keeping  this  symmetry  in  mind  one  finds  that  the  values  of 
the  roll  moments  due  to  bow  seal  (FKBS) ,  stern  seal  (FKSS), 
propulsion   (FKP)   and  plenum  pressure  (ABPB*PHI*Z)  are  also 
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symmetric  with  respect  to  zero  moments.  For  propulsion  and 
plenum  pressure  the  sign  of  the  moments  are  identical  to 
that  of  the  roll  angle  while  for  bow  and  stern  seal  mcments 
the  signs  are  reverse.  From  the  results  given  in  TABLE  VII 
with  positive  steady  state  roll  angles  Figure  4  has  been 
drawn  where  the  dotted  line  represents  the  craft's  initial 
state,  the  dashed  line  represents  the  craft's  steady  state 
considering  the  transom  deadrise  angle  and  the  solid  line 
represents  the  craft's  steady  state  (smaller  roll  angle) 
considering  the  center  of  gravity  deadrise  angle  for  the 
vertical  force  calculation.  Same  markings  apply  to  the 
moment  vectors  shown  from  which  their  change  in  magnitude 
for  both  deadrise  angles  can  be  seen. 


*     SW 


Figure    4    -       ROLL    MOMENTS    IN    A    PORT    TURN 
(ROLL    MOMENTS    AND    ANGLE    NOT    TO    SCALE) 


Legend:         AED   Aerodynamic 
ES      Bowseal 
SS      Sternseal 
RUD   Rudder 
Deadrise    anaia: 


CG   Center    of   gravity 
PP    Plenum    pressure 
SW    Side walls 
P      Propulsion 

-  -   Transom 

Center    of    gravity 
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Sines  the  simulation  run  with  the  center  of  gravity 
deadrise  angle  being  used  for  the  vertical  force  computation 
gave  results  that  were  more  favorable  (better  damped  time 
histories  and  steady  state  values  closer  to  those  measured 
experimentally)  this  change  in  vertical  force  computation 
was  implemented  into  the  XR-3  Loads  and  Motion  Program. 
Starting  with  the  deadrise  angle  given  in  Table  7  for  draft 
up  to  seven  inches,  again  straight  line  approximations  for 
the  effective  deadrise  angle  at  the  center  of  gravity 
location  (station  14)  for  larger  drafts  have  been  performed 
resulting  in  an  almost  linear  relationship.  This  part  of 
subroutine  SIDEWALL  is  now  valid  for  any  draft  and  contains 
the  following  statements  which  have  been  used  in  this  form 
for  all  simulation  runs  to  follow  in  this  study: 

DDRANG=0.0 

IF  (DS.GT. 0.5833)  DDRA NG=  (DS-0 . 53 33)  *0 . 0629 

DRANG=1 .021+DDRANG-PM1*PHI 

CTNDF=COTAN  (DRANG) 

FZHOLD  (J) =FZH  (J) 

FZHDRP  (J)=PM1*FYH  (J)  *CTNDR*  PROMO  1 

FZH  (J)  =FZH  (J)  +FZHDRP  (J) 


FK=  (FZH  (2)  -FZH  (1)  )  *  YS  W+FKD-FY*ZS 

It  should  be  noted  that  vertical  projected  forces 
from  both  sidewalls  are  used  in  computing  the  roll  moments 
in  the  XR-3  Loads  and  Motions  Program.  Section  III .  E  of  this 
thesis  will  be  concerned  with  the  effect  of  changing  the 
simulation  program  to  consider  the  deadrise  projected  force 
only  from  the  outward  sidewall  in  a  turn  maneuver  as  it  is 
discussed  in  Sect.  III. 3. 2  . 
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C.   CROSS-FLOW  DRAG 


In  the  preceding  section  the  lateral  force  FYD  which  is 
due  to  cross-flow  drag  has  been  found  to  be  rather  large 
compared  to  FYHf  thereby  changing  the  direction  of  the 
resultant  force  significantly.  FYD  is  computed  in 
subroutine  SIDEWALL  by  summing  the  contributions  DF  of  all 
stations : 

DF(I,J)  =  -HRH0*CDSW*VREL*A3S(VREL) *DSWAV  (I,J) 

where    OF  (I, J)  is  the  lateral  drag  contributed  by 

the  j-th  element  of  each  (i-th)  sidewall 
HRHO  is  RHO/2. 
7REL  is  relative  velocity 
DELX  is  the  incremental  distance  used  for  the 

wetted  draft  calculation  along  each  sidewall 
DSWAV  is  the  average  wetted  draft  over  one 

incremental  distance 
CDSW  is  the  cross-flow  drag  coefficient  for  one 

sidewall  . 

The  cross-flow  drag  coefficient  being  used  in  the  present 
XR-3  Loads  and  Motions  Program  is  the  value  corresponding  to 
that  for  flew  normal  to  a  long  flat  plate  being  CDSW=1.28  . 
Ref.  3  points  out  that  an  investigation  of  cross-flow  drag 
forces  and  coefficients  (before  the  vertical  added  mass  A33 
had  been  reformulated  as  discussed  in  Sect.  III.A.1)  showed 
better  agreement  between  tneoretical  and  experimental 
results  if  CDSW  =  2.t),  representing  the  normal  drag 
coefficient  for  a  sharp  flat  plate,  had  been  used  instead. 
After  the  reformulation  of  A33  the  agreement  between 
modified   theory   and   experiment  had  very  much  improved,  so 
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that  the  former  drag  coefficient  could  be  used  again.   Since 

in   this   thesis   work  the  change  in  A33s  did  not  effect  the 

roll  behavior  significantly  (TABLE  IV)  ,    some  turn   maneuvers 

at  20  kn  and  15  degrees  rudder  angle  have  been  simulated  for 

various  drag  coefficients  in  order  to  investigate  which  CDSW 

would   better   match  simulated  and  experimental  steady  state 

roll  angle.   From  TABLE  VII  a  decrease   in   roll   angle   and 

roll   moment  is  desireable  and  therefore  the  cross-flow  drag 

force  FYD  should  be  decreased  via  reducing  CDSW.  The   steady 

state   values  are  shown  in  TABLE  VIII  with  the  lateral  added 

mass  coefficient  being  unchanged  C  =0.4. 

h 

TABLE_VIII 

Steady  state  conditions _f or_various_CDS j 

20  kn  turn,  15  deq  port  rudder  angle 

CDSW  1 . 28        1.16         1.0         0.7 


Pitch  angle 

0.43 

0.43 

0.44 

0.46 

Roll  ar 

igle 

1.55 

1.48 

1  .38 

1.14 

Draft 

5.73 

5.80 

5.83 

5.91 

Speed 

19.39 

19.36 

19.31 

19.19 

FKBS 

-343.0 

-330.5 

-308. 1 

-252.8 

FKSS 

-   2.7 

-   2.6 

-   2.4 

-   2.0 

FKSW 

207.8 

175.0 

127.8 

40.0 

FKRUD 

-  82.2 

-  47.9 

-   2.3 

80.4 

FKAED 

-  56.6 

-  58.7 

-61.6 

-  67.6 

FKP 

0.  14 

0.14 

0.  13 

0.  1  1 

AEPB*PHI*Z 

276.4 

264.6 

246.5 

201.9 

FYH  (P, 

S) 

-4  3/-125 

-45/-125 

-48/-124 

-57/-121 

FYD 

-654.9 

-616.9 

-561 . 0 

-436.0 

Comparison  with  Ref.  6  (Table  VI)  shows  close   agreement 
in  roll  angle  (1.36°)   for  the  run  with  CDSW=1.0  .  But  these 
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trial  runs  have  been  executed  with  arbitrary  smaller  drag 
coefficients  because  an  increase  in  drag  coefficient  (e.g. 
CDSW=2.0  for  sharp  flat  plate  as  in  Ref.  3)  would  have 
enlarged  the  disagreement  between  simulated  and  experimental 
data.  Also,  from  plot  29  for  CDSW=1.28  and  plot  31  for 
CDSW=1.16  it  is  seen  that  a  smaller  drag  coefficient  results 
in  more  damping  in  the  roll  response,  in  this  case  six 
versus  nine  cycles  of  distinguishable  oscillations.  Since 
the  choice  of  CDSW=1.0  cannot  be  declared  to  be  appropriate 
for  the  actual  shape  of  the  XR-3  sidewalls  it  will  not  be 
considered  any  further  in  this  study.  CDSW=1.16  generated  a 
4  %  change  (about  15  %  is  desired)  toward  the  experimentally 
measured  roll  angle  and  its  effect  will  be  investigated 
together  with  the  present  cross-flow  drag  coefficient 
CDSW=1.23   in  the  studies  to  follow. 

Shown  below  in  Figure  5  is  the  dependence  of  steady 
state  roll  angle  in  turns  generated  by  a  15  degrees  rudder 
angle  to  port  at  20  kn  on  the  cross-flow  drag  coefficients 
used  in  the  simulation  runs  (TABLE  VIII)  . 


Roll  angle 


1.5" 

1.4- 

1.3- 
1.2-- 
1.1.. 


1.0- 


0.7 


0.8 


0.9 


1.0 


1.1 


1.2 


-^CDSW 


1.3 


Figure  5  -   S.3.  ROLL  ANGLE  VERSUS  CDSW 
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D.   THRUST  MAPPING 


In  turn  maneuvers  as  they  have  been  discussed  in  the 
preceding  sections  another  important  point  to  be  considered 
is  the  change  in  thrust  on  both  engines.  Since  the  engines 
mounted  on  the  XR- 3  provide  constant  power  for  a  given 
throttle  setting,  the  actual  thrust  delivered  during  a  turn 
maneuver  may  vary  even  though  the  throttle  remains  fixed. 
In  all  simulation  runs  shown  in  this  thesis  until  now 
constant  thrust  has  been  used.  In  Ref.  6  experimentally 
measured  reduction  in  thrust  is  given  for  turn  maneuvers  at 
20  kn  and  various  rudder  angles.  From  this  and  Ref.  7 
providing  more  detailed  information  the  thrust-rudder 
dependence  was  as  listed  below. 

TABLJLIX 
Reduction  in  .thrust  fpr  various  rudder  angles 

in.  turn  maneuvers,  at  2  0  _kn 

Rudder  angle      total  reduction  in  thrust  {Ref.  6) 
(deg)  {%) 

5  -  2.8 

9  -  3.3 

12  -4.3 

TA3LE  IX  has  been  used  for  thrust  mapping  (Block  16  of 
the  data  input)  in  the  simulation  runs.  Due  to 
unavailability  of  more  precise  data  it  has  been  assumed  that 
the  total  reduction  in  thrust  is  contributed  from  both 
engines  in  egual  amounts.  3ut  this  may  not  be  guite  true, 
since  due  to  the  craft's  roll  angle  a  difference  in 
immersion  depth  for  both  engines   could   result   in   unegual 


changes  in  thrust.   TABLE  X  shows  the  steady  state  values  of 

these    simulations    considering   two    cross-flow    drag 

coefficients,    deadrise    force   contributions   from   both 

sidewalls  computed  at  the  center  of  gravity  as  well  as   both 

lateral  added  mass  coefficients  C  =0.4(0.8). 

h 

TABLE_X 
Steady  state  ..conditions  using  thrust  mapping 
in  20  kn  -turns_at_yarious_rudder_anqles 


Rudder 

angle 

(deg) 

CDSW=1 . 

,28 

5 

9 

12 

CDSff=1. 

,  16 

5 

9 

12 

Roll  angle  (deg) 
simulation 


0.30(0.25)  0.09 
0.71(0.61)  0.28 
0.99(0.39)    0.58 


0.28(0.23)  0.09 
0.67(0.57)  0.28 
0.94(0.34)    0.58 


Speed  (kn) 
simulation   Bef.  6 


19.3(19.2)  19.3 
18.9(19.0)  19.5 
18.5(13.6)    19.2 


19.2(19.2)  19.3 
18.9(19.0)  19.5 
18.4(18.6)    19.2 


Comparing  the   steady   state   roll   angles   obtained   if 

thrust   mapping   is   used   with  those  listed  in  TABLE  VI  (no 

thrust  mapping,  deadrise  force  computation  at  the  center   of 

gravity) ,  a  change  in  roll  angle  by  6  to  1 4  %  for  C  =0.4  and 

h 

by  14  to  20  %   for  C  =0.8  toward  the  measured   angle   can   be 

h 

observed.    The   steady   state   roll   angles  found   in   the 

simulation  runs  still  differ  from  the  roll  angles  given  in 
Ref.  6  by  a  factor  of  about  two  and  three  but  can  be 
expected  to  be  a  little  closer  to  these  if  proper  thrust 
mapping,  i.e.  different  maps  for  port  and  starboard  engine, 
is  used.   By  this  also  the   steady   state   velocities   which 
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differ   by   about   2.5  %   can   be  expected  to  come  in  tetter 
agreement,  i.e.  to  increase  for  smaller  roll  angles. 

Hegarding  the  roll  angle  time  responses  for  simulation 
runs  using  both  lateral  added  mass  coefficients  (0.4  and 
0.8)  and  12  degrees  rudder  angle  the  following 
characteristics  have  been  found:  the  runs  using  thrust 
mapping  (Table  X,  plot  33(35))  show  almost  the  same 
transient  behavior  as  do  the  runs  without  thrust  mapping 
(Table  VI,  plot  21  (23))  with  about  ten  cycles  of  oscillation 
before  the  transients  die  out.  But  the  runs  using  thrust 
mapping  need  quite  a  long  time  (t>45  sec)  before  they  reach 
the  steady  state  roll  angle. 

The  effect  of  using  a  larger  lateral  added  mass 
coefficient  (0.8)  can  be  observed  in  plot  35.  Comparing 
this  with  plot  33 ,  the  roll  angle  response  in  plot  35  shows 
a  smaller  peak  roll  angle  (1.25°  versus  1.47°)  and  the 
transients  have  smaller  amplitudes.  But  both  dia  out  after 
ten  cycles.  This  effect  has  already  been  noted  in  Section 
III.B.4. 

Investigating  the  effect  of  changing  the  coefficient 
CDSW=1.28  to  1.16  by  comparison  of  the  roll  responses  shown 
in  plots  33  and  37  or  35  and  39,  a  significant  damping 
effect  as  has  been  found  previously  in  Section  III.C  can  not 
be  observed  here. 
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E.   DEADRISE  FORCE  OF  OUTWARD  SIDEWALL 


This  part  of  the  investigation  on  how  to  effect  the  roll 
behavior  of  the  testcraft  simulation  is  again  concerned  with 
the  deadrise  angle. 

In  the  present  XR-3  Loads  and  Motions  Program  deadrise 
forces  for  both  sidewalis  are  computed  with  the  deadrise 
angle  appropriately  corrected  by  the  roll  angle  for  port  and 
starboard  sidewalis.  The  idea  behind  this,  e.g.  in  a  port 
turn,  is  to  consider  a  pressure  buildup  (large  FIH(S))  on 
the  outboard  side  of  the  starboard  sidewall  and  a  redaction 
in  pressure  on  the  outboard  side  of  the  port  sidewall  (small 
FYH(P)).  Thereby  an  upward  vertical  force  FZHDRP  (with 
negative  sign)  is  obtained  for  starboard  side,  while  a 
downward  vertical  force  (positive  sign)  is  computed  for  port 
side.  Both  vertical  forces  are  then  added  to  the  buoyancy 
force  of  their  side.  The  forces  obtained  by  this  approach 
are  shown  in  Figure  6. 

In  Ref.  4,  however,  a  different  philosophy  in  regarding 
the  acting  forces  is  presented.  As  discussed  in 
Section  III.E.2  an  approximate  check  on  the  craft's  roll 
stability  can  be  made  if  the  SES  geometry  and  the  vertical 
location  of  its  center  of  gravity  are  known  (Fig.  1).  In  a 
port  turn  the  principal  force  effecting  the  roll  behavior 
arises  from  the  outboard  side  of  the  starboard  sidewall 
(relative  high  pressure  on  the  structure  there  due  to  wave 
buildup)  while  there  is  only  a  small  force  (due  to  small 
pressure  change)  on  the  inboard  side  of  the  port  sidewall 
contributing  a  roll  moment  only  if  there  exists  a  roll 
angle.  Now  both  vertical  forces  are  directed  upward.  This 
approach  and  the  resulting  forces  are  shown  in  Figure  7. 
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FYH(P)« 


FZHDRP(P) 


Figure  6  -   ACTING  DEADRISE  FORCES  (TWO  SIDEWALLS) 


FYH(P) 


Figure  7  -   ACTING  DEADRISE  FORCES  (ONE  SIDEWALL) 

Note  :  Since  the  vertical  force  contributed  by  the  inboard 
sidewall  is  rather  small,  as  sho*n  above,  it  can  be 
neglected  in  the  roll  stability  check  (Sect.  III.B.2). 
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To  investigate  the  effect  of  this  interpretation  of  the 
acting  forces  two  statements  controlling  the  computation  of 
deadrise  forces  via  the  sign  of  the  rudder  angle  RUDANG 
(port  rudder  is  positive)  and  the  sidewall  under 
investigation  (PM1=-1  for  port  side)  have  been  added  to  the 
computation  procedure  for  the  vertical  projected  force  : 

DDEANG=0.0 

IF     (DS.GT. 0.5833)     DDRANG=  (DS-0 . 58 33)  *0 . 0629 

DRANG=1.021+DDRANG-PM1*PHI 

CTNDR=COTAN  (DRANG) 

RUDSIG=SIGN (1 . , RUDANG) 

IF  (RDDSIG.NE. PM1)     CTNDR=PH1 *T AN  (PHI) 

FZHOLD(J)  =FZH  (J) 

FZHDRP  (J)  =PH1  *FYH  (J)  *CT NDR*PROMO 1 

FZH  (J)  =FZH  (J)  +FZHDRP  (J) 

The  steady  state  values  for  simulated  turn  maneuvers  using 
thrust  mapping  and  center  of  gravity  deadrise  angle  are 
shown  below  considering  two  cross-flow  drag  coefficients  and 

both  lateral  added  mass  coefficients  C  =0.4(0.8). 

h 

TABLJLXI 
Steady  state  conditions 
in ,20  kn  turns_ at  various  rudder  angles 


Rudder 

ang 

le 

(deg) 

CDSW=1. 

,28 

5 

9 

12 

CDSH=1 , 

.16 

5 

9 

12 

Roll  angle  (deg) 
simulation   Ref.  6 


0.43  (0.  52) 

0.09 

0.90(1 .00) 

0.28 

1.22(1.  34) 

0.61 

0.41(0.51)  0.09 
0.87(0.98)  0.23 
1. 17(1.30)    0.61 


Speed  (kn) 
simulation   Ref.  6 


19.4 (19. 5) 

19.3 

19. 1 (19. 4) 

19.5 

18.7 (19.0) 

19.  2 

19.3(19.5)  19.8 
19. 1  (19. 3)  19.5 
18.6 (19.0)    19.2 
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Comparing  the  steady  state  values  for  simulated  runs  in 
tables  X  and  XI  ,  the  effect  of  using  only  the  outward 
sidewall  for  the  deadrise  force  computation  can  be  observed 
to  give  less  favorable  (little  larger)  steady  state  roll 
angles  but  mere  favorable  (larger)  steady  state  velocities. 
This  is  exactly  what  has  been  expected  since,  as  shown  in 
Figures  6  and  7,  if  deadrise  force  contributions  from  both 
sidewalls  are  used,  the  deadrise  force  from  the  port 
sidewall  will  reinforce  the  one  from  the  starboard  sidewall 
toward  a  less  outward  roll  moment  (Fig.  6)  while  applying 
the  philosophy  depicted  in  Fig.  7  the  vertical  force  from 
the   port  sidewall  counteracts  the  starboard  deadrise  force. 

Regarding  the  roll  angle  responses  shown  in   Appendix   A 

for  12  degrees  rudder  angle  with  CD5tf=1.28  (plots  41  and  43) 

and  CDSW=1.16  (plots  45  and  47)  and  comparing  these  with  the 

corresponding   ones   from   the  previous  section  (plot  33,  35 

and  37,  39)  the  change  in   the   vertical   force   computation 

considering   the   outward  sidewall  only  is  seen  to  result  in 

about  8  versus  10  cycles  of  transient  oscillations  with  less 

amplitude.   Using   C  =0.8   (plots   43   and   47)   has  a  quite 

h 

significant  damping  effect  and  makes  the  roll  angle  response 

more  resemble  a  step,  which  is  in  accordance  with  the 
experience  gained  by  Sef.  7.  Also,  comparing  these  plots 
with  plots  35  and  39,  the  negative  roll  angle  occurring  at 
the  moment  when  the  rudder  has  been  introduced  does  not  show 
up  any  more. 
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F.   VERTICAL  LOCATION  OF  CENTER  OF  GRAVITY 


In  Section  III.B.2  it  was  discussed  how  to  check  for  the 
testcraft's  roll  bahavior  in  a  turn  provided  the  5SS 
geometry  and  the  location  of  the  center  of  gravity  are 
known.  This  section  will  be  concerned  with  the  effect  of 
different  vertical  locations  of  the  center  of  gravity  on  the 
simulated  XR-3  roll  behavior 

From  experimental  measurement  the  location  of  the  center 
of  gravity  of  the  XR-3  has  been  determined  to  be  10.05  ft 
forward  of  the  transom  and  Leo  and  Boncal  established  the 
vertical  location  at  2.54  ft  above  the  keel  on  the 
longitudinal  center  line  of  the  craft.  As  modifications  to 
the  testcraft  were  introduced  (e.g.  engines  and  seals)  the 
horizontal  location  of  the  center  of  gravity  was  again 
established  experimentally.  No  such  changes  were  made  on  the 
vertical  location.  In  order  to  investigate  the  sensitivity 
of  the  steady  state  roll  angle  to  vertical  locations  of  the 
center  of  gravity,  reductions  in  height  (in  accordance  with 
the  stability  check  geometry  depicted  in  Figure  1)  in 
increments  of  0.  1  ft  were  entered  into  the  Loads  and  Motion 
Program. 

Simulation  runs  with  constant  thrust  and  initial  speed 
of  20  kn  introducing  a  15  degrees  port  rudder  angle, 
deadrise  force  contributions  from  both  sidewalls  computed  at 
the  center  of  gravity,  the  sidewall  drag  coefficient 
CDSW=1.28  and  the  lateral  added  mass  coefficient  being  0.4 
were  executed  and  the  results  are  shown  in  TaDle  XII.  From 
the  listed  data  it  follows  that  with  a  lower  vertical 
location  of  the  center  of  gravity  a  tendency  toward  an 
inward  roll  angle  (here   actually   a   smaller   outward   roll 


45 


angle)  can  be  achieved.  As  shown,  a  decrease  in  vertical 
distance  by  0.1  ft  resulted  in  about  0.1  degree  decrease  for 
both  peak  and  steady  state  roll  angle  (almost  linearly 
related) .  From  the  corresponding  plots  as  indicated  in 
Table  XII  a  roll  damping  effect  due  to  the  relocation  of  the 
center  of  gravity  can  hardly  be  recognized.  The  number  of 
cycles  of  transient  oscillations  and  the  steady  state  speed 
were  about  the  same  for  all  runs  listed  in  Table  XII. 

TABLE_XII 
Roll  conditions  in  20  kn  turn  with_constant  _t^r^§^ 
and  1-5  degrees _pprt  rudder  _anqle 

vertical  location  of  CG      experimental 
2.54   2.44   2.34   2.24     testcraft  data 

steady  state 

roll  angle     1.55   1.43   1.32   1.22         1.36 
peak 

roll  angle     1.67   1.53   1.46   1.35 
peak 

roll  rate      2.13   2.07   1.98   1.91 
cycles  of 

transients     8      8      8      8 
refer  to 

plots         13-14  49-50  51-52  53-54 
steady  state 

speed         19.4   19.3   19.3   19.3  18.7 

Table  XIII  shows  the  effect  of  changing  the  vertical 
location  of  the  center  of  gravity  for  simulation  runs  at 
20  kn  with  12  degrees  port  rudder  angle  and  thrust  mapping 
applied,  deadrise  force  computation  for  both  sidewalls  at 
the  center  of  gravity  and  the  coefficients  being  0.4  for  the 
lateral  added  mass  and  CDSW=1.28  for  the  sidewall  drag. 
These   conditions   were  similar   to   those  used  in  Table  X. 
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From  Table  XIII  the  effect  of  reducing  the  vertical  location 
of  the  center  of  gravity  in  0.1  ft  increments  over  a 
reasonable  range  can  be  observed  to  also  result  in  linear 
reductions  in  maximum  and  steady  state  roll  angle  for  the 
case  if  thrust  mapping  is  used.  Nonlinear  reductions  appear 
for  the  peak  roll  rate.  The  number  of  cycles  during  the 
transient  period  (ten)  is  again  not  effected  by  the  vertical 
relocation  of  the  center  of  gravity.  Regarding  the 
corresponding  plots  a  very  slight  damping  effect  can  be 
observed  for  reduced  vertical  locations.  The  results  given 
in  Table  XIII  are  graphically  represented  in  Figure  8  . 


Ti3LE_XIII 
Roll  conditions  in  20  kn ^urn^with^thrust  mapping 
and  12  degre^s^port  rudder  ..angle 

vertical  location  of  CG      experimental 
2.54   2.44   2.34   2.24     testcraft  data 

steady  state 

roll  angle     0.99   0.91   0.34   0.77  0.53 

peak 

roll  angle     1.47   1.36   1.25   1.16 
peak 

roll  rate     2.14   1.92   1.74   1.58 
cycles  of 

transients    10     10     10     10 
refer  to 

plots         33-34  55-56  57-58  59-60 
steady  state 

speed         13.5   18.4   18.4   13.4  19.2 
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Figure  3  -   ROLL  CONDITIONS  VERSUS 
VERTICAL  LOCATION  OF  CENTER  OF  GRAVITY  (TABLE  XIII) 
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G.   ROLL  EAMFING  DDE  TO  VERTICAL  WAVE  GENERATION 


One  part  of  subroutine  SIDEWALL  is  concerned  with  the 
roll  damping  effect  which  is  due  to  vertical  wave  generation 
during  a  roll  motion.  The  computation  of  the  roll  damping 
term  was  developed  by  Oceanics  Inc.  and  added  to  the 
program  by  Menzel  [  Ref .  2].  The  principles  used  in  the 
development  cf  this  addition  were  not  provided  by  Oceanics 
Inc.,  however,  the  apparent  idea  is  to  reduce  the  computed 
roll  moment  FK  by  some  value  which  is  dependent  on  the 
craft's  roll  rate  and  its  draft,  thereby  on  the  vertical 
added  mass  A33.  During  the  study  of  the  effect  of  these 
additional  terms  a  discrepancy  was  noted  in  the  dimensions 
used  in  the  calculations. 

The  roll  moment  FK  is  calculated  by  the  equation  given 
in  SDWL1950  and  considering  the  expressions  leading  to  this 
statement  it  is  found  that 

dim  [FK]  =  lb  .  ft 

which      is      the      correct      dimension      for    a    moment   vector.    Now 
considering   the   dimension   of    the    roll    damping   terra 

dim   [PR0M02*¥SW*YStf*BC2*P/PI  ] 

which    is   subtracted    from   FK    (SDWL2200)     it    is    found 

PRCM02   and    PI    are    dimensionless 

dim   [YSW  ]    =    ft 

dim  [P]   =  rad/sec 

dim  [BC2  ]  =  dim  [  AC2] 

=  dim  [FC2  *  length] 

=  dim  [F2  *  cos(x)  *  length] 

=  dim  [  A33  *  length] 
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=  dim  [  RHO  *  B  *  B  *  length] 

=  (lb  .  sees  /ft*)  .  ft  .  ft  .  ft  . 

So  for  the  roll  damping  term 
dim  [YSW  *  YSW  *  P  *  BC2] 

=  ftz  .  (rad/sec)  .  (lb  .  secVft*)  .  ft* 
=  lb  .  ft  .  sec 

which  is  not  the  correct  dimension  for  a  moment  expression. 
Assuming  that  only  the  given  terms  may  be  involved  in  the 
roll  damping  calculation,  it  is  possible  to  come  up  with  the 
proper  dimension  if  the  roll  rate  P  were  squared  but  keeping 
its  sign,  thus  getting 

roll  damping  term  =  PROH02*YSW*YSW *BC2*P*ABS (P) /PI  . 

Also  the  vertical  added  mass  should  be  associated  with  some 
•velocity  squared'  -  expression  which  is  absent  in  the 
original  version  of  this  part  of  the  program  but  could  be 
generated  by  squaring  the  roll  rate  P.  The  effect  of  the 
modified  expression  will  be  an  increase  in  the  damping 
effect  for  roll  rates  P>1.0  which  exist  only  in  the  initial 
phase  of  the  turn  maneuver  (maximum  roll  rate  in  all 
simulation  runs  was  about  2.0)  .  But  for  most  of  the  run 
length  P  is  less  than  unity  and  the  supposed  modification 
will  decrease  the  damping  effect  especially  when  approaching 
steady  state  (P=0.0)  . 

The  PR0MG2  -  term  provides  a  means  to  arbitrarily  set  a 
roll  damping  factor  which,  as  stated  in  Reference  2,  has 
been  found  from  experiments  to  be  16.0  .  Until  now 
PROMO  2=1.0  has  been  used  in  this  thesis  work.  So  the  effect 
of  changing  it  to  the  experimental  value  is  investigated 
next.  In  the  runs  to  be  discussed  turn  maneuvers  at  20  kn 
with  a  12  degrees  port  rudder  angle  have  been  simulated  with 
deadrise  forces  being  computed  at  the  center  of  gravity, 
using   the   experimental  damping  factor  PROMO2=16.0   and  the 
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coefficients    CDSW=1.28    for    the  sidewall    drag   and   C    =0.4      for 

h 

the  lateral  added  mass. 


1.  Thrust  is  held  constant  with  deadrise  force  contributions 
from  both  sidewalls 

Comparing  the  obtained  plots  61-62  with  plots  21-22 
(PROMO2-1.0)  only  a  slight  damping  effect  on  roll  angle 
response  can  be  observed  (peak  roll  angle  1.45°  versus 
1.47°),  but  the  number  of  transient  oscillations  decreases 
from  eleven  to  nine.  The  peak  roll  rate  decreased  from  2.14 
to  2.10  for  PROMO2=16.0  and  the  damping  effect  is  more 
pronounced  in  the  roll  rate  time  history. 

2.  Thrust   is   mapped  with  deadrise  force  contributions  from 
both  sidewalls 

Comparing  the  plots  63-64  with  the  corresponding  ones 
(33-34)  for  PROM02=1.0  ,  again  only  a  slight  damping  effect 
can  be  noticed.  The  peak  roll  angle  is  again  reduced  from 
1.47  to  1.45°  and  the  number  of  cycles  of  transient 
oscillations  is  uneffected. 

If  the  damping  term  is  changed  to  include  the  P*ABS(P) 
expression  (plots  65-66)  there  result  eleven  versus  nine 
cycles  of  transient  oscillations  with  little  higher 
amplitudes  (peak  roll  angle  is  again  1.47°  as  it  was  for 
PROMO2=1.0)  .  The  roll  rate  time  history  also  shows  little 
larger  amplitudes  (peak  roll  rate  is  2.14  versus  2.10)  . 

3.  Thrust  is  mapped  with  deadrise   force   contribution   from 

the    outward    sidewall   only   and   lateral   added   mass 

coefficient  c  =0.8 
h 

Comparing   these   plots   (67-68)   for  PROi«I02  =  16 . 0  with  plots 

43-44  where  PRCMO2=1.0   again  only  a  slight   damping   effect 

in   roll   angle   response   can  be  noticed  with  the  number  of 
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transient  oscillations  being  uneffected.  Using  the  P*ABS(P) 
term  the  slight  damping  effect  shown  in  plot  67  is  cancelled 
as  shown  in  plot  69. 

The  above  observations  show  that  the  experimentally 
determined  roll  damping  factor  ?ROMO2=16.0  has  only  a 
negligible  damping  effect.  Considering  the  roll  damping  term 
again  as  it  was  added  to  the  present  XR-3  Loads  and  Motions 
Program  in  1974,  there  was  a  1/PI  -  factor  which  actually 
cancelled  the  PI  -  factor  contained  in  the  former  expression 
for  the  vertical  added  mass  A33S  (Section  III. A).  Since  the 
PI  -  factor  already  has  been  eliminated  by  correcting  the 
A33S  -  expression,  there  is  no  further- need  to  cancel  it  in 
the  damping  term  by  the  1/PI  -  term  which  therefore  also 
could  be  eliminated.  This  change  is  expected  to  cause  an 
increase  in  sidewall  roll  damping  by  a  factor  of  PI  and  has 
been  investigated  by  an  additional  run. 

4.  Thrust  is  mapped  with  deadrise   force   contribution   from 

the  outward  sidewall  only,  roil  damping  term  being 

16.0*YSW*YSW*BC2*P 

and  lateral  added  mass  coefficient  C  =0.8 

h 

From  plots  71-72  again  only  a  slight  damping   effect   versus 

plots  67-68  can  be  noticed. 

The  investigation  described  in  this  section  shows  that 
the  roll  damping  due  to  vertical  wave  generation  as  included 
in  subroutine  SIDEWALL  is  not  very  effective  even  if  the 
factor  PR0M02  in  the  damping  term  is  increased  from  16.0  to 
50.0  .  There  also  is  a  lack  in  matching  dimensions 
(statement  SDWL2200)  . 
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IV.   PROPULSION  AND  RUDDER  SUBROUTINES 


In  all  simulation  runs  executed  during  this  study  the 
XR-3  Loads  and  Motions  Program  as  given  in  Ref.  2  has  been 
used  except  for  the  modifications  mentioned  in  this  work. 
This  simulation  program  contains  Forbes'  version  [Ref.  9] 
of  ' SUBROUTINE  PROP1  which  included  some  revisions  of  the 
original  version  given  by  Leo  and  3oncal  [Ref.  8].  Since  in 
TABLE  VII  (Section  III.B.4)  the  roll  moment  contributed  by 
the  propulsion  system  (FKP)  is  rather  small  (0.1  to  0.2)  an 
additional  run  using  Leo  and  Boncal's  version  has  been 
executed  under  the  conditions  as  stated  for  TABLE  VII 
(center  of  gravity  deadrise  angle).  The  steady  state  values 
are  shown  below  (in  parantheses  are  the  values  obtained 
using  Forbes1  version)  : 


pitch  angle    0.44(0.45)     draft      5.73(  5.81) 
roll  angle      1.67(1.53)     speed     19.39(19.37) 


FYSW 

-909.  4(- 

■821.9) 

FYAED 

-  21.3  (- 

■  20.0) 

FYP 

97.  0( 

0.0) 

FKBS 

-348.  0(- 

■336.4) 

FKSW 

269.  9  ( 

212.  1) 

FKAED 

-  63. 5(- 

■  56.0) 

ABPB*PHI*Z 

299. 5( 

273.5) 

FYRUD   -  54 .7  (   32.7) 
3*V*AM   888.4  (  809.2) 


FKSS  -  3.0(-  2.7) 
FKRUD  150. 3  (-  89.9) 
FKP     -305. 2(   0.14) 


FYD  -729. 2(-651. 4)    FYH  (P/S)   -43/-1 37 (-45/-1 26) 
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Regarding  these  values  drastic  differences  are  found  for 
FYP  and  FYRUD  as  well  as  for  FKP  and  FKRUD,  though  the 
steady  state  values  for  pitch  and  roll  angle,  draft  and 
speed  are  close  to  those  obtained  using  Forbes'  version. 
The  pitch  and  roll  angle  responses  are  given  in  plots  25-28. 
Comparing  these  with  the  corresponding  ones  using  Forbes' 
version  (plots  13-16)  it  is  seen  that  Leo  and  Boncal's 
version  generates  more  oscillatory  roll  and  pitch  angle 
responses.  Therefore  a  review  of  subroutine  PROP  especially 
the  computation  of  the  respective  forces  and  moments  has 
been  carried  out  and  are  discussed  below. 

The  review  resulted  in  force  and  moment  eguations  which 
are  identical  to  those  given  by  Leo  and  Boncal  [Ref.  8]. 
Another  point  that  supports  the  supposition  to  use  this 
version  of  subroutine  PROP  is  the  fact  that  it  provides  a 
reasonable  roll  moment  due  to  the  propulsion  system  (Forbes1 
version  provides  almost  zero  roll  moment) .  In  order  to  tsst 
the  effect  of  changing  subroutine  PROP  to  its  original  form, 
simulation  runs  1.  and  2.  from  the  preceding  section  have 
been  repeated  and  their  time  histories  for  the  roll  behavior 
are  shown  in  plots  73-74  and  75-76,  respectively.  Comparing 
these  with  the  corresponding  ones  using  Forbes'  version  of 
subroutine  PROP,  plots  61-62  and  63-64,  it  is  observed  that 
the  steady  state  and  peak  roll  angles  are  larger  by  about 
7  %  and  the  number  of  transient  oscillations  is  uneffected 
by  this  change. 

Next  a  review  of  subroutine  RUDDER  has  been  carried  out 
in  order  to  investigate  the  cause  of  the  change  in  sign  for 
FYRUD  and  FKRUD  depending  on  the  propulsion  subroutine  used 
(Forbes1  or  Leo  and  Boncal's  version)  although  identical 
RUDDER-subroutines  have  been  used  in  both  cases.  The 
important  eguations  to  be  considered  are  : 
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DSR=Z+ZS-XE*THETA 

ENDFAC=1.    +    DSR/(DSR  +  RSPAN) 

VH=V    +    XR*R    -    ZR*P 

EFFANG=RUDANG  -  ENDFAC*VH/U 

FY=RHO*U*U*RAREA*ENDFAC*RCLB*EFFANG 

• 

FK=  -ZR*FY 

The  above  equations  compate 

-  draft  at  the  rudder  location 

-  endfactor   depending   on   the   rudder   length   below   the 
craft's  keel 

-  sideslip    velocity 

-  effective   rudder   angle  depending  on  the  craft's  forward 
and  sideslip  velocities 

-  lateral   force   depending   on   rudder   area  and  effective 
rudder  angle 

-  roll  moment  due  to  rudder  . 

The  sign  of  the  roll  moment  FK  is  dependent  on  the  sign 
of  the  lateral  force  FY  which  depends  on  the  sign  of  the 
effective  rudder  angle  EFFANG.  This  in  turn  depends  of 
course  on  the  introduced  rudder  angle,  but  the  magnitude  of 
the  second  term  of  this  equation  (sideslip  angle)  is 
responsible  for  the  actual  sign  of  EFFANG.  VH/U  represents 
the  tangent  of  the  sideslip  angle  which  -  applying  small 
angle  approximation  -  is  subtracted  from  the  introduced 
rudder  angle  (see  Figure  9).  Since  ENDFAC=1.73  and  0  was 
the  same  for  both  runs,  a  difference  in  VH  which  was  larger 
in  Leo  and  Boncal's  version  by  about  50  %  caused  the 
negative  sign  (versus  a  positive  sign  using  Forbes'  version) 
for  the  effective  rudder  angle  EFFANG  and  thereby  FY.  So 
the  only  way  to  introduce  a  change  toward  an  inward  roll 
behavior  here  (FKRUD  was  positive)  is  to  first  check  the 
input  value  for  RSPAN  and  second  to  eventually  modify  the 
craft's  engines  to  have  larger   immersion   into   the   water. 
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With  regard  to  the  first  item,  if  RSPAN  would  increase, 
ENDFAC  would  decrease  as  would  the  magnitudes  of  EFFANG 
(negative)  ,  FY  (negative)  and  FK  (positive)  . 


SSA 


ESSA 


>VH 


EFFANG 


SSA  =  VH/D 

ESSA=ENDFAC*VH/U 

RUDANG 

EFFANG 


k-^  RUDANG 


ESSA 


sideslip  angle 
effective  sideslip  angle 
rudder  angle 
effective  rudder  angle 


Figure    9  -   EFFECTIVE  RUDDER  ANGLE  IN  A  TURN  MANEUVER 


Considering  next  the  fact  that  the  engines  of  the  XR-3 
had  been  replaced  in  1976  a  recent  measurement  of  the  rudder 
and  propeller  location  yielded  the  following  data  (original 
input  values  shown  in  parantheses)  : 


XPO  =  -1 . 50 
ZPO  =  -1 .333 
XRO  =  -1. 083 
RSPAN=  1.333 
RAREA=  1.42 


-1.  275) 

-0.604) 

-  1.  12  5) 

1.21  ) 

0.68  ) 
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Two  of  the  runs  discussed  in  the  preceding  section  have 
been  repeated  using  the  original  roll  damping  terra  and  the 
new   propulsion   and   rudder    input    parameters. 

1.  Thrust   is    mapped    with    deadrise    force      contributions      from 

both    sidewalls    with   C    =0.4    . 

h 

Comparing  plots  77-78  with  the  corresponding   ones   for   the 

former   rudder   and   propulsion  parameters  (plots  63-64)  the 

expected  effect  due  to  a  decrease  in  FKRUD  can   be   observed 

resulting   in   smaller   peak  roll  angle  (1.34°  versus  1.45°) 

and  a  significantly  decreased  steady  state  roll  angle  (0.65° 

versus  0.99°  at  t=50  sec) .   The  transient  period  of  the  roll 

angle  response  in  plot  77  shows  the  same  number  of  transient 

oscillations   as   plot   63   but   the   response  did  not  reach 

steady  state  even  after   45  sec  .    The   steady   state   roll 

angle  could  be  estimated  to  be  about  0.3  to  0.4  degrees. 

2.  Same   conditions  as  in  the  previous  run  but  with  deadrise 

force   contribution   from   the   outward  sidewall  only  and 

C  =0.8  . 
h 

Regarding  plots  79-80  versus  plots  71-72  a  similar  effect  as 

mentioned  before  is  observed.  Now  the   roll   angle   response 

shows   a   more   oscillatory  transient  period.   The  peak  roll 

angle  reduced  from  1.46  to  1.36  degrees  while  the  roll  angle 

at  t=50  sec  changed  from  1.34  to  0.97  degrees.   For  this  run 

the  steady  state  roll  angle  could  be  estimated  to   be   about 

0.8   degrees   since   the   roll   angle  response  again  did  not 

reach  steady  state  for  the  run  using  the  recently  determined 

rudder    and   propulsion   parameters.    The   initial   small 

negative  roll  angle   which   appeared   in   the   previous   run 

(plot  77)    at    the    moment   when   the   rudder   motion   is 

introduced,  does  not  show  up  any   more.    The   overall   roll 

response   shown   in   plot  79   does   not   fall  off  as  steeply 

(compared  to  plot  77)  if,  in  a  turn  maneuver,  deadrise  force 

contribution  from  the  outward  sidewall  only  is  considered. 
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From  these  runs  it  is  obvious  that  the  rudder  and 
propulsion  parameters  of  the  new  engines  result  in  a  more 
favorable  XR-3  roll  behavior  with  a  significantly  less 
steady  state  outward  roll  angle.  At  the  time  of  this 
sensitivity  study  Reference  7  pointed  out  that  there 
actually  is  no  roll  angle  indicated  by  the  measuring  devices 
in  a  port  turn  at  20  knots  and  12  or  15  degrees  rudder 
angle.  Keeping  the  device  sensitivity  of  ±0.5  degrees  in 
mind,  it  can  be  concluded  that  the  steady  state  roll  angle 
provided  by  the  present  XR-3  Loads  and  Motions  Program 
including  the  before  mentioned  program  modifications  by  the 
author,  is  quite  well  in  agreement  with  the  actual  craft 
data.  Further  efforts  should  be  undertaken  to  obtain  a 
better  damping  effect  during  the  transient  period. 

Remark  : 

Although  the  new  engines  extend  about  0.5  ft  deeper  into  the 
water  than  the  old  engines,  the  new  input  value  for  RSPAN  is 
close  to  the  former  value.  In  Reference  1  the  denominator  of 
the  second  term  in  the  equation  for  ENDFAC  is  defined  to  be 
the  distance  from  the  water  surface  to  the  bottom  tip  of  the 
rudder.  The  equation  for  ENDFAC  used  in  the  present  program 
accounts  for  changes  in  draft  DSR  in  both  the  numerator  and 
the  denominator.  The  RSPAN-term  should  then  be  defined  to 
be  the  distance  from  the  craft's  keel  to  the  bottom  tip  of 
the  rudder.  The  author  suspects  that  formerly  there  has  been 
chosen  too  large  a  value  for  RSPAN  in  the  case  of  the  old 
engines  (RSPAN=1.21  ft,  probably  from  water  surface  to 
bottom  tip  of  the  rudder) .  Since  the  bottom  tip  of  the 
rudder  and  the  center  of  the  propeller  were  located  at 
almost  the  same  distance  from  the  craft's  keel  the 
magnitudes  of  ZPO  and  RSPAN  should  be  about  the  same  which 
is  not  true  for  the  former  chosen  values.  From  these 
observations  the  author  concludes  that  until  now  a  too 
optimistic  roll  behavior  resulted  due  to  a  too  large  value 
for  RSPAN. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 


The  preceding  sections  of  this  chapter  reflect  the 
results  of  an  investigation  on  how  the  XR-3  roll  behavior  in 
turn  maneuvers  at  20  kn  simulated  with  the  Loads  and  Motions 
Program  is  affected  by  certain  changes  in  parameters  and 
forces.  The  results  of  this  sensitivity  study  are 
summarized  as  follows  : 

*  Added  mass  effect 

A  reformulation  of  the  vertical  added  mass  computation  which 
reduced  the  A33s-value  by  a  factor  of  1/PI  had  only  a  little 
effect  in  changing  steady  state  draft,  pitch  angle  and 
required  thrust  for  constant  speed  in  straight  runs.  In 
turn  maneuvers  this  change  resulted  in  less  damped  responses 
for  pitch  and  roll  angle  and  affected  the  steady  state 
values  for  pitch  angle  by  +8.3  %  ,  roll  angle  by  -1.0  %  and 
draft  by  +  1.4  %  .  Considering  the  roll  behavior,  the  effect 
of  reducing  the  A33s-value  by  a  factor  of  1/PI  in  the 
vertical  added  mass  computation  was  negligible. 

A  change  in  the  lateral  added  mass  coefficient  C  =0.1  to 

h 

0.8    had   the  effect   of   reducing   the   amplitudes  of  the 

transient  oscillations  and  the  steady  state   roll   angle   by 

0.1   degree  toward  the  experimentally  measured  value.   Since 

in  a  turn  maneuver  with  15  degrees  rudder  angle  an   unstable 

craft   behavior   showed  up  which  had  not  been  experienced  in 

practice,  the  value  of  C  =0.8  is  not  realistic   if   deadrise 

h 

force  contributions  from  both  sidewalls  are  considered.  If 
deadrise    force   contribution    from    the   outward      sidewall      only 
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is  used  in  the  simulation  program  the  increase  in  lateral 
added  mass  coefficient  causes  an  increase  in  the  steady 
state  roll  angle  by  about  0.1  degree.  But  since  now  the 
roll  angle  response  is  quite  significantly  damped  and 
resembles  a  step  as  it  does  in  practice,  the  choice  of  the 
larger  coefficient  (0.8)  seems  to  be  appropriate  for  the 
case  of  the  deadrise  force  contributed  from  the  outward 
sidewall  only. 

*  Deadrise  angle  effect 

The  deadrise  angles  over  the  length  of  the  curved  XR-3 
sidewalls  are  not  uniform.  Changing  the  deadrise  angle  which 
is  used  in  the  deadrise  force  computation  from  the  transom 
(78.8°)  to  the  center  of  gravity  (58.5°)  resulted  in  more 
damped  responses  for  pitch  and  roll  motion  in  turn 
maneuvers.  For  all  chosen  rudder  angles  generating  port 
turns  at  20  kn  the  steady  state  roll  angle  changed  favorably 
by  23  to  38  %  toward  the  corresponding  roll  angle  measured 
experimentally  and  documented  in  Ref.  6. 

*  Cross-flow  drag  coefficient 

Changing  the  cross-flow  drag  coefficient  (CDSW=1.28  for  a 
long  flat  plate)  to  lower  values,  e.g.  CDSW=1.00  ,  the 
testcraft's  simulated  roll  behavior  could  be  forced  to  match 
the  measured  steady  state  values.  But  since  this 
coefficient  has  been  chosen  arbitrarily  and  cannot  be 
declared  to  be  appropriate  for  the  actual  XR-3  sidewalls 
(CDSW=1.0  means  that  there  is  no  effective  drag  coefficient) 
this  CDSW-value  has  not  been  considered  any  longer.  Another 
suggested  cross-flow  drag  coefficient  CDSW=1.16  [Ref.  3] 
resulted  in  about  5  %  reduction  in  steady  state  roll  angle 
and  a  shorter  transient  period  (about  20  %)  with  smaller 
amplitudes   of   oscillation    than   did   the    runs    using    CDSW=1.28. 


*  Thrust  mapping 

I 
Performing   thrust   mapping  in  the  simulation  runs  using  the 

data  provided  in  Ref.  6  and  assuming  that  both  engines   lose 

equal   amounts   of   thrust   in  turn  maneuvers,  the  simulated 

steady  state  roll  angles  could  be  shown  to  agree  better  with 

the  measured  angles  by  about  another  10  %  but  the  roll  angle 

response  did  not  quite  reach   steady   state   even   after   45 

seconds. 

*  Deadrise  force  from  outward  sidewall  only 

The   simulation   runs   executed   with   the   deadrise   force 

computed  for  the  outward  sidewall  only  and  thrust  mapping  as 

done   before,   showed   less   agreement  for  both  C  =  0.4  (0.3) 

h 

with  the  corresponding  measured   steady   state   roll   angles 

than   did   the   runs   considering   deadrise  forces  from  both 

sidewalls.   But  for   C  =0.8   the   roll   angle   response   was 

h 

significantly  damped  and  looked  more  like  a  step,  as  was 
experienced  in  practice. 

*  Vertical  location  of  the  center  of  gravity 

The  given  value  for  the  vertical  location  of  the  center  of 
gravity  (2.54  ft  above  the  keel)  might  not  be  quite  correct 
due  to  the  difficulties  in  experimentally  determining  it. 
The  significant  effect  of  the  vertical  center  of  gravity 
location  on  the  XR-3  roll  behavior  is  obvious  from  Figures  1 
and  8  and  the  discussion  in  Section  III.?.  Since  it  controls 
exceedingly  the  X3-3  roll  behavior  it  is  important  that 
future  investigators  obtain  an  accurate  measurement  of  the 
C.G.  location  before  taking  testcraft  data. 


*  Subroutine  PROP 

A  review  of  the  propulsion  subroutine  revealed  that  the 
force  and  moment  calculations  as  given  by  Leo  and  Boncal 
(Ref.  8)  should  be  used  instead  of  the  version  stated  by 
Forbes  (Ref.  9).  This  change  results  in  a  7  %  larger  outward 
steady  state  roll  angle  but  it  provides  a  reasonable  roll 
moment  due  to  the  propulsion  system. 

*  Roll  damping  due  to  vertical  wave  generation 

The  roll  damping  calculation  included  in  subroutine  SIDEWALL 
has  only  a  negligible  effect  on  the  craft's  net  roll 
behavior  over  the  range  of  the  damping  factor 
1 . 0<PROM02<50.0  .  Since  there  also  is  a  discrepancy  in  the 
dimension  of  the  roll  damping  term  the  source  of  this  part 
of  the  program  should  be  investigated. 

*  Suggested  areas  for  further  investigation 

During  the  studies  undertaken  with  the  XR-3  Loads  and  Motion 

Program   the  author  found  the  following  points  either  in  the 

input  data  from  the  testcraft  or  in  the  program  that  require 
further  investigation  : 

■  Rudder/thrust  angle 
As  learned  from  Ref.  7  the  rudder/thrust  motion  is  directly 
introduced  to  the  port  engine  while  the  starboard  engine 
follows  this  via  a  metal  bar  connecting  both  engines.  The 
starboard  rudder  and  thrust  vector  may  be  off  by  2  or  3 
degrees  to  either  side,  depending  on  the  direction  of  the 
rudder/tnrust  action.  Thereby  the  roll  behavior  will 
certainly  bs  effected.  Since  rudder/thrust  mapping  is 
available  in  the  Loads  and  Motion  Program,  it  is  important 
that  accurate  data  be  collected  for  each  engine  during  a 
test  run. 
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■  XR-3  weight 

During  each  test  run  there  is  a  possibility  that  the  weight 
is  not  exact  since,  as  learned  from  Ref.  7,  the  testcraft 
takes  on  water  in  the  seals  and  thereby  its  weight  increases 
causing  larger  draft  which  also  effects  the  roll  behavior. 
This  additional  weight  can  be  substantial  and  attempts 
should  be  made  to  measure  or  estimate  the  additional  weight 
to  be  used  in  the  Loads  and  Motion  Program. 

■  Cross-flow  drag  coefficient 

Further  efforts  should  be  undertaken  to  establish  the  proper 
value  of  the  cross- flow  drag  coefficient  since  CDStf=1.28 
might  not  be  quite  appropriate  for  the  actual  shape  of  the 
XR-3    sidewalls. 

■  Sideslip    velocity 

Those  parts  of  the  Loads  and  Motions  Program  where  terms 
contributing  to  the  lateral  velocity  V  are  computed  should 
be  reviewed  since  a  large  value  of  V  may  result  in  a  change 
in  the  effective  rudder  angle  from  positive  to  negative. 

In  this  sensitivity  study  it  has  been  shown  that  the 
most  significant  effect  on  the  simulated  XR-3  roll  behavior 
are  contributed  by 

*  deadrise  angle 

*  vertical  location  of  the  center  of  gravity 

*  rudder  and  propeller  location  . 

Less  significant  effects  are  contributed  by 

*  thrust  mapping 

*  cross-flow  drag  coefficient  . 


63 


For  further  studies  concerned  with  the  XR-3  Loads  and 
Motions  Program  for  improved  roll  behavior  representation  it 
is  recommended  to  use 

-  the   reformulated  equation  for  the  vertical  added  mass 
A33 

-  deadrise  angle  at  the  center  of  gravity  location 

-  deadrise  force  contribution  from  the  outward   sidewall 
only 

-  subroutine  PROP  as  given  by   Leo   and   Boncal   and   in 
Appendix  C  of  'this  thesis 

-  thrust  and  rudder  mapping. 
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APPENDIX  A 


PLOTS 


In  the  table  given  below  are  listed  all  plots  with  the 
respective  parameters  used  in  the  corresponding  simulation 
run.  The  table  contains  the  following  abbreviations  : 
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pitch    rate 
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PITCH     DQTF 


D 


u  - 1 0 

r, 

NG 

RQ 

•  >  F  R .  j  i 

1  .1 

T  T  ! 

,  r~ 
1 1— 

i  ; 

PLOT     12 
for    applied    parameters    see    first    page    of    this    appendix 


80 


b 


!  I 


K-SCRLF--d-OOE+Oi   UHIT5    INCH. 
'-SCRlE   5-OOE-Oi  UNITS    IMC1-. 

RGR0G5    ■ 
PLOT    IS 


TURN   20   KN 
ROll   RNGl- 


RUDM=15 

r~  Pi  p-  :   i  p- 


T1 


PLOT    13 
for    applied    parameters    see    first    page    of   this   appendix 


W\A 


mv^ 


v^ 


121 


K-3CRLE-i.0GE+0i   UNITS    INCH, 
---5rqLc-5.03E-0:  units  inch. 

RGR0D5    .    TURN   20   KN 
PLOT    IS   ROLL   RRTE 


RUDM=15 

VERSUS   time 


PLOT     14 
for    applied    parameters    see    first    page    of    this   appendix 


032 


12: 


33L 


v*-5C;RLF.=4.00E-0i 

RGR0D5    p    T 

plot  is  ph 


UNITS 
UNITS 

JRN 
~H   F 


INCH. 
INCH. 


on 


m 


w 


p. 


r-- 


,s 


!  irnOi  If       T  T 


r- 


PLOT    15 
for    applied    parameters    see    first    page    of    this   appendix 


33 


mkhtAWA 


V  v  v  V  w 


A  A  A 


K-3CRLE  --1.  OOE+Oi   UNITS    INCH, 
r-5r-RLF--2.Gnr.-n2  UNITS    INCH, 

RGR0Q5    ,    TURN   20   KN    .    RUDM-iR 
PLOT    IS        PITCH   RRTF  UER5U5 

PLOT    16 
for    applied    parameters    see   first    page    of    this   appendix 


T  IMF 


84 


A/\AAAAA/WW\ 


mm 


K-SCflLF.^1.  OOF.+Oi  UNITS  INCH. 
-'-5rGu---5.00P-'Ji  UNITS  INCH. 

RGR0E5    ,    turn  20  KM 


PLOT    13        ROLL   RN 


r~ 


u 


Riin=i  n 


VJ 


NO 


DP 


ERSUS   TIME 


PLOT    17 
for    applied    parameters    see    first    page    of    this    appendix 


aaamaamMa 


WvYt 


K-SCRLE>1.00E+0i  UHI^s    !MCH> 
v-5CRl5'--5.  JGE-Gi  UNITS    IHCH* 

RGR0E5    ,    TURN 


D 


lGT 


t  -: 


n 


ROLL 


20 
RRTE- 


KN 


RUO-iO    =    HO   RQ 
UERSUS    TIME 


PLOT    18 
for    applied   parameters    see    first    page    of   this    appendix 


IIAMa 


X-Sr.RLF^  l-OCE+Oi   JNITS    INCH 
V-SCRLF~5.QGE-0i   UNIT5    INCH 


r  n  r 


1  <J 


PL3T 


h  ! 

13 


TUR 

ROLL 


20 


K 


hGl:: 


I    irp; 

Or  K 


w><     O     O 


TTMF: 


PLOT     19 
for    applied   parameters    see    first    page    of   this   appendix 


87 


3J- 


PLOT 


l  j 


RHLL    RHGLc 


UERSUS   TIME 


K-5CRLE-J...  OOr-O:  UNITS  INCH. 
V-5CRLE--2.  GOF-G1  UNITS  INCH, 

PLOT    21 
for    applied    parameters    see    first    page    of   this    appendix 


89 


I  ft 


k-sr.rlf»1.0ce+qi  jkits  inch. 
■scrl;  -o-oce-qi  units  inch. 

.      T!  IRK 


,fjRT 


ni 

t- 


R 

L3T    15   ROLL   RRTE 


n 


JiJUJ 


Q      T  T  ML" 


PLOT    2  2 
for    applied   parameters    see    first    page    of   this   appendix 


90 


llAA/v 


V 


i:j 


SrflLE- 


LOOF.+Oi  UNITS    INCH 
2.00F.-01   UNITS    INCH 


RGRT63 
PLOT    15 


TURN   20   KN» 
ROLL   RNGLE 


,  irn 


T    T  k,j  ~ 


PLOT    2  3 
for    applied    parameters    see    first    page    of   this   appendix 


91 


\/Vw 


Tr5*~r 


K-3CRLF.-4-  GOE+Oi  UNITS  INCH 


-SCRlE 

RGRT63 


. OOE-Ol  JMIT5    INCH, 

-    TURH   20   KN 


n 


l_U 


HT 


R( 


H 


RRTE 


UER3U3 


T  M 


PLOT  24 
for  applied  parameters  see  first  page  of  this  appendix 


92 


f<     '  ■  AAA 


;   r   •' 


>0i   JNITS    INCH- 


■Cm     -  .-.     "Gj.      U 


-1-  :      i  i—    l     ! 

li.   w 

!/! 

i — i  , — * 

- 

1        ' 

J    ! 

:l 


(—■ 


r 


r  ^ 


PLOT    2  5 


for    applied   parameters    see    first    page    of   this   appendix 


93 


II  ft  H , 

JJldJl/l  L  A  A  A  ,A  ;■■  ^i^w^  r  ,-, »  ^~,  ^ 


-  -\  i  u_i  i 


^  °  p  .j.  n 


s  i  ki 


!  ;  r 


D  r;  P f 

Di    DT     T 


j  -  OCr.  -J.i   JNj  !  b    lK;,.h 
Tj  !Pki    9fl 


KN» 


LU 


h  U  l  L 


n  n  t  c 
n.r_i  !  l 


J     h :  J 


T  ki  C 


PLOT    26 
for    applied    parameters    see    first    page    of    this   appendix 


94 


\ 


'  y  If !  1 1)  v  * y ' 


K-5tqi_c  i.^Or-j;.  UNITS  INCH 
('-SCRLr-  1  ncr-Gi  UNITS  inch 


-.  p.  n  D  .  3 


m  i 


20   KH 


PLOT    15   PITCH   HNGl-: 


ia 


RQ 


;  i  p  n  ~- 
[y  L  h  -_j 


Tit 


kir 


PLOT    2  7 
for    applied   parameters    see   first    page    of   this    append 


IX 


95 


r1-: 


Cs     ! 


!       I 

il 


,j  !  i 


i 


i  i 


!  ! 


l!  J 
i  f 


— 


il    ! 


i  !    I 


il 


I 

!! 


I       I:      , 

i    I,    ! 

I 


irrtt- 


U;  I  i  i 


li 


1  A 


I 

:l       I 


I  i 


I   !!   I  I    I 


li    I 


:i 


rhT 


_U- 


I  i 


I 


li    v 


M 


:  T   T 


j     i 


1"  I 


i   i  i  * : 


H'^r       no     I    ki  ;t.'.;      Tki^i 


PLOT    2  8 
for    applied    parameters    see    first    page    of   this   appendix 


96 


L" 


A/yvw- 


n  i 


'KITS    TkPH 


INCH 


r- 


K-3rRL_E,-.L  QGH>0 

v-5r.RLF  -3-0C?-0i   JN1T 

RGRTS1    ,    TURN   20   KK 
PLOT    IS        RGll   RNGLE 


ijFR 


v^  O 


THF 


PLOT    2  9 
for    applied    parameters    see   first    page    of    this   appendix 


97 


D 


Uu 


Pr^- 


35: 


r  1/ 


K-j 


?CRLF 


1«  GCE+OI  UNITS  INCH. 
■01  UNITS  INCH, 


■J  ■  LL 


RCRT3i    ,    TURN   2D   KH 
^lQT    IS   ROLL   RRTF 


TT 


PLOT    30 
for    applied    parameters    see    first    page    of    this   appendix 


98 


~1 


uL 


A/vw 


K-3CQLE-:L0CF>0i  JHIT5    INCH. 

RGRT82    ,    Tj 


l_-  ^.L'ur-ui  UNITS    INCH, 


D 


RN   20    KH 
lQi    IS        RQ^l   RNGLE 


;  I  p  n  c  i  I  c  T  T  ki  P 


PLOT  31 
for  applied  parameters  see  first  page  of  this  appendix 


99 


(\    A 


Ut 


:y  v  ^ 


X-SCRLF. ---I.  OOF.+OI  JNITS    INCH- 


srn 


LT 


IN  IT: 


nrpT  n  ") 


DL 


in: 
j    Hull    hh ! t 


TURN 


^0   KM. 


UER3US   TIMF 


PLOT  32 
for  applied  parameters  see  first  page  of  this  appendix 


100 


A  A 


nc 


PLOT    is       roll   RNGlE 

K-5CRLr>i.0DF+0i  UNITS    INCH- 
V-SCRLF^2«00F.-Qi   JNIT5    INCH- 

PLOT    33 
for    applied    parameters    see    first    page    of   this   appendix 


UcRjUj   TIME 


101 


K-SCfiLE  ■■•i-COF+Oi  UNITS    INCH. 

^-SC^l.-  -5.  DOF-Oi  UNITS    n-.'rH, 

RGRTKi    -    TURN   20   KH. 
DlQT    IS   ROll    RRTE 


UERMU5   TIME 


PLOT    34 
for    applied    parameters    see    first    page    of    this    appendix 


102 


■A 
o 


a 

-t 
n 


121 


K-SCRLrd-OOE+Oi  UNITS  INCH. 
'-'  -SCRLC  -'2 .  OCF-0 1  JH 1 T5  !  NCH . 

RGRTK2    , 
PLOT    IS 


TURN   20   K 
ROlL   n.HGLE 


;  i  p  r  c  i  I  q  t  T  ^  r2 


PLOT    35 
for    applied    parameters    see    first    page    of    this   appendix 


103 


fflll  V 


^VW- 


1   V 


'-'-5rRL-- 


?0F+01  UNITS  INCH 
i  INCH 


nr!-  n  i  ,  'K'  t 

UUu.   U-l   UH  x 


RGRr*2    ,    TURN   20   KN, 
-LOT    IS   ROll   RRTE 


n 


l  r~"  P;  c-     i  r~-       t  r 


PLOT    36 
for    applied    parameters    see    first    page    of    this    appendix 


104 


,;-5rPLr  l-OCF+Oi  UNITS  INCH. 
v-Sr.qL~  2.0CE-Q1  UNITS  INCH. 

RG 

pl  a 


RGRT:<3    ?    TURN 


I J 


i  r 


l_  n  u  ^  sj 


roll  rngle;         I 

PLOT    3  7 
for    applied    parameters    see    first    page    of   this   appendix 


TIM 


105 


lux  i' 


¥ 


7\rv< 


i 


K-5CRLr..---1.00E."+01  UNITS    INCH, 


UNITE 


"  -J-  .HLl.-  -J  •  UUl  -(Jo. 

RGRTK3    ?    TURn   ^ 
PLGT    13   RQLL    RRTE 


TkT.hi 


k 


i    I  i —  H  r-  I    IP        T   T  I  ,-r  r" 

iJEFbuj      I!  • 


PLOT    3  8 
for    applied    parameters    see    first    page    of    this   appendix 


106 


< 

-5r 

?i~ 

i 

,DCF.+Qi 

Ulli    1  -J 

INCH.. 

,1 

-Sr 

JL~ 

n 

,  OGE- 

-Q  ! 

Ui  1  i     1    ^J 

INCH. 

RG 

RT 

x' 
r  i 

4      > 

T 

1  Ipkl 

o  r  t    i 

£-L>       : .    : 

N> 

— ^ 

■  • 

'lj: 

T 

I 

-•> 

R 

Oll 

RNGL 

i — 

irn 


■J  U  J 


T  M  F 


PLOT  39 
for  applied  parameters  see  first  page  of  this  appendix 


107 


K-5CPLF  -1 .  OOF- 01  JN 1 T3  1 NCH . 
-'-SrcL-  -5.0C- '-01  UNITS  INCH. 

RGRTK4    ,    T'jRN   20   KN 
DLQT    T3   ROLL   RRTE 


UER3U3   TIME 


PLOT    40 
for    applied   parameters    see    first    page    of    this    appendix 


108 


t: 


T  ' 


PLOT    13        ROLL    RN 


ji_r 


UER5US 


T   l/if 


kiU- 


^-Sr^Lr.--i^GG--Gi  UNITS  INCH: 
V-SCRLF^-OOF-ai  UNITS  INCH. 

PLOT  4  1 
for  applied  parameters  see  first  page  of  this  appendix 


109 


I 


l/l  A  A 


i^UU 


!  fctt  1 1  i  v  J 


TV v 


1-OGE+Oi  UNITS  INCH., 

00F.-01   UNITS    INCH, 
t:  mki     on     [^ 


•'-5rPL-  3 


RGRT 

PLOT 


T      -~ 

i.J 


RH 


L 


R 


u 

i— - 

r 


Jrh 


!  r 


h  .1  """ 


JU-J 


PLOT    42 
for    applied    parameters    see   first    page    of    this    appendix 


110 


.  rv-"W\/v 


- 


HT     "  ~- 


ROLL   R 


—  -.• 


_  - 


•j1     i —    ■  _  - -  _.    ■_  1  _  i  -j     uiij 


-j  n  c  _  r  ■ ,  -  -  -      F  K  -r  M 

•_    _i  -  -        U-i-         ■_  j.         \J  -         vl .        ' 


PLOT    43 
for    applied    parameters    see    first    page    of   zhis   appendix 


111 


lift 


gywwwfcr? 


'vV    -/^l^^i-^^— '^W , 


i:. 


K-3CRLE-1. OOF+01  JHIT3  INCH. 
{■-SrfiLF^.OOF-Oi  UNITS  INCH, 

RGRT12    ,    TURN   2G   V 
PlQT    13   RQlL   RRTF 


UERSUS    TIME 


PLOT  4U 
for  applied  parameters  see  first  page  of  this  appendi: 


112 


V 


y  \l  IAAa 


nc2 


DLQT    IS 


RGuL   RNGL 


UERSU3   TIM 


n-5rRL=  i^DC~+0:  UNITS  INCH, 
r-5nRLE-2,GGF.-0i  UNITS  INCH. 

PLOT    45 
for    applied    parameters    see    first    page    of   this   appendix 


113 


; 


11 A 


\   A   A 


pnM 


rr^= — srr 


<-scrle;-i-Gl?e+ui 


i  1 


^ 


SCRLE-5 

RGRT1 
PlGT    [ 


III1! 


J 


NITS 
11  UNITS 


INCH. 
INCH, 


T 


U- 


RH   20   K 


RG 


L 


RRTE 


UERS 


i  r 


U 


T  T  ML 
!    1  !  i  l 


PLOT    46 
for   applied    parameters    see    first    page    of   this   appendix 


114 


!j2 


12. 


PLQT    15        ROLL    RHGLE 

K-5C,RLE==1.00E+0i  UNITS  INCH. 
^-5CRLE.-=2.0GE-01  UNITS  INCH. 


!  |  p  n  C I  I  C      T  T  U  p 


PLOT    4  7 
for    applied   parameters    see    first    page    of   this    appendix 


115 


W\/V\AA^? 


■■/v^.rv*^Tr 


V" 


i:. 


v 


X-SHRLF.^l.  OGF-^OL   UNITS    ^NCH. 
y-SrRLF.--5-GCF.-0i   JNITS    INCH, 

RGRT14    ,    TURN   20   KM- 
PL3T    13   ROlL   RRTE 


UERSU5 


TM 


PLOT    48 
for    applied    parameters    see    first    page    of   this   appendix 


116 


33* 


32E, 


•:« 


dlQT   IS       ROlL  rngle 


VERSUS    TIME 


X-SCHlE; inOOE+OI  UNITS  INCH. 
V-5CCILE^2.00E-01  UNITS  INCH, 

PLOT    4  9 
for    applied    parameters    see    first    page    of   this   appendix 


117 


1  i  M  ,  ,  ,'  > 


b* 


i  :  r-i_ 


K-3GRLE--  1.  DOE-Oi  UNITS 
^-3CnLE=  5.  Du^-UA  UNITS  INCH- 


RUN0Q1 


D 


□T 


i     i —  i 


rURN   20   KN 

RRTF 


L/C 


r-  p-  n- ,   I  r- 
H   ~i!   J   1 


T  THF 


1 1  i 


PLOT    50 
for    applied    parameters    see    first    page    of   this    appendix 


118 


dlQT 


lc- 


3j2 


3  J- 


R3ll 


^GlE 


UERSUS   TIME 


fr-SCRlX-I- OQE+Qi   UNITS    IH 
^-SCRlE--  Z.OCE-QI  UNITS    Ih 


PLOT    51 
for    applied    parameters    see    first    page    of    this   appendix 


119 


A/V^1  -  -g 


3J., 


K~SCRLE»L.0GE+Q1  UNITS  INCH- 
V-SCRLE^OOE-Qi  UNITS  INCH. 


RUN0D2    .    TURN   20 
DL0T    is  roll   rrte 


KK 


UERSUS   TI 


1 F 


PLOT    52 
for    applied    parameters    see   first    page    of   this   appendix 


120 


X-5CRLE-1.DQ"E 
V-SCRLE----2 

RUN0D3 


-01 


INCH. 


nntr     'i  i 


u-l    iJNITS 


UNITS 

INCH, 

TURN  20  KH 
RQll   ANGLE 


UERSUS   TI 


PLOT    53 
for    applied    parameters    see    first    page    of   this   appendix 


121 


■fljj 


3  J 


K-5CRLE--  1-  G0C.J-01  JH:ts  IHC-H- 
?-5CRLF  5.00E-01  UH!ts  ikCH, 

RUN0D3    -    TURN   2u   kk 
PLOT    IS   RGl.l    RRTE 


FRr 


1 


t  i^  r-" 


PLOT    54 
for    applied    parameters    see    first    page    of   this    appendix 


122 


l  VvVv 


PLOT 


:i 


ROLL 


1         |l         t    -^J     I L-.    . 


iJFR 


K-5Ccu^  1   ?2[r+0i   uNITS    INCH- 
^-5rRLF.--2.QCF.-01  u^ITS    INCH. 

PLOT  55 
for  applied  parameters  see  first  page  of  this  appendix 


123 


'  IL-  T  TQ      T  i  '■  r-i; 


K-5CRLE-.L  dOE-^Ol 

/-3CRLF:  5.  OuE-Oi  UNITS    IN! 

RGRG- 1     -    tjjrkj   20 
PLOT    IS   RQl.l   RRTE 


KM 


n;  in   1  '  1 

H  U  LJ  "  1  L1 
:  r  h  "j 


t  T  Mr 


PLOT  56 
for  applied  parameters  see  first  page  of  this  appendix 


124 


oj: 


3J.J 


033 


i:« 


X-SCRLE=i.OOE+Oi  UNITS  INCH. 
V-SCRLE--2.00E-01  UNITS  INCH. 


D 


RGRQF4 

LQT  is 


TURN  20   K 
ROLL   ANGLE 


RUO=iO    ,    NO   RQ 


UER 


C 


5 


TI 


f 


PLOT    57 
for    applied    parameters    see    first    page    of    this   appendix 


125 


3C4 


K-SCFLE-I.OOE+ai  UNITS  INCH, 
y-3CRLF.--5.00E-0i  UHITS  INCH. 


RRR0P4 

PLOT    15 


TURN   20   KN    .    RUD-10    ,    HQ   RQ 


ROLL   RRTE 


ri  ,J  u  _-• 


i    TTMF 


PLOT  58 
for  applied  parameters  see  first  page  of  this  appendi: 


126 


V  V  \A 


I  I        \J  »  •  ■    |  L.  .-  ■     J.  «    v.1  |J  l  .  X  (J  J.        Uill     I   J  -  I  T  V (         1 


D 


r^i 


ifiE\q 


PLQT 


r  r 


Ti  !RN 

*  :ULl 


\/ 


Li  NJ  i     i     — 


H  U  u ; 


r 


r  t  k,j  r 


PLOT    59 

for    applied    parameters    see    first    page    of   this    appendix 


127 


r-^-£ 


1/ 


-frr- 


vvvy 


RLh 


01  UNITS  INCH, 
DOE -01  JN.IT5  INCH. 


RHRilFq 


j 


Tl  iQkj 


KN 


n  i  n  t     t  q    n  n  i 
r  L  U  !      I  J    H  J  L 


L_ 


hh  !  r 


r>  ;    i 


!UDM=J 


_y 


JU 


T  I M f 


PLOT    6  0 
for    applied    parameters    see    first    page    of    this    appendix 


123 


PLOT    I 


r 


i—.  n ,    ;        o  k  ;  r 


•n1   UNITS   INCH. 


JFRS'JS    T"Ih' 


V  L_ 


.SfPL!    2,0Cr,-01   JNITS    tNCH 


for  applied  parameters  see 


PLOT  61 

first  page  of  this  appendix 


129 


A  A 


\  V 


K-SCRLE ■■  1.  OOF  -Oi  UNITS  INCH. 

-j-5C^lF-5.0CF-01  UNITS  INCH. 

RGR0E'.6    .    TL!RN    2Q    KH    -    RUD-lrj 

RRTE 


PLOT    15   RGll 


I  !  P  R  q !  : ~ , 


N 


RR 


i  _J 


/]L 


PLOT    6  2 
for    applied    parameters    see    first    page    of    this   appendix 


130 


PLOT    IS 


i — 


RQu-L   RNGll 


VLri  Jo  J       !   1 !  I 


^-Sr.qLE-i.GGF+O.i.  UNITS    I 
v-Sr;RLB-:2.00F.-QI   UNITS    t 


NGH 


1U1    I  ■> 


PLOT    63 
for    applied   parameters    see   first    page    of   this   appendix 


131 


K-SCRLE:  i-OGE+0. 
Y-SCRLE^-3.  DGE  -0. 


u 


k:TTS 


r.H, 


!HITS    INCH. 

PLOT    13   RGl.l   RRTE 


dtp  n  ^~  q 

n  'o  n  L  L  y 


TURN 


D 


i  c: 

l 


L_. 


T  f  M c 


PLOT    54 
for    applied    parameters    see    first    page    of   this    appendix 


132 


I  (I  • 


i:. 


PLOT    15        ROuL   RNGlE 


UERSU5   TIME 


K-SCRLi::--l-QOE+01   JNITS    rNCu. 
v'-5CRLF.-2^CF.-Qi   lJMt3    INCH, 


PLOT    6  5 
for    applied    parameters    see    first    page    of   this   appendix 


133 


7°Cto-? 


K-5t',RLE^l.GDE+0i  UNITS    INCH. 

•--5r-L?-5-0ur-Oi   JNIT5    ?HCH. 


RGRDL9    ,    TURN   2G 
PLOT    13   RGl.l   RRTr 


K 


m  !  n  -  i 


RliR 


•^ 


R3US 


T  T  k/f 
I   i    i 


PLOT    66 
for    applied    parameters    see    first    page    of    this    appendix 


134 


PLOT    15 


ROLL  rngle: 


vthbuo     '  liit 


f i      J '  ,  ■    1—  .        J.»  uti.  '  u.1.      U !  i  i    I  -J  l  "kw    i  ' 


T!SCH 


PLOT    67 
for    applied   parameters    see    first    page    of   this    appendix 


135 


v:' 


u 


x 


■^■/'■//'/■M^v^ 


K-5r:RLF    I.:Cr-Gi  UNITS    INCH, 
^-5CRLr   5.0CE-0I  J!-:T5    tHC-H. 

^GRC:"n    ■    TURN   20   KM 


n; 


::l.HT    IS    R( 


n 


L_  L_ 


1 .    I 


r 


i  i  n  -  i  r- 

ULj  —  _l  J 

:  !rnr-vi  i 

I  -  !      '  l      J  :    J 


K 


i  n     n  n 


|     1^/1  ! 
II' 


PLOT    68 
for    applied    parameters    see    first    page    of    this    appendix 
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<A/V\ 


plot  is      rulL  rn: 


r~" 


:  I  rr  n  c» , 


JL. 

K-SHRLF    1.  nOF-Hi  lIHITS    INCH. 
^-Sf^LF   2.GCE-01  UNITS    INCH. 

PLOT    69 
for    applied    parameters    see    first    page    of    this   appendix 


3    T  [  k 
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UT 


i^n/vVl/V^vv- 


•■•V  '-'.V.'--  ■-»  ■-■  ' 


^-srqi--- 

.  I 

lJKITS 

INCH 

■  0  ■ 

n  °  r     o  ' 

UNITS 

ThC^ 

PfjDDFR  T!  1RN    ?H    KK  CMP-1 


Pi   PIT     IS    PH!   i      RDTP" 


v  L.  n  vJ  U 


T  T  Ivi  p" 
1    I  I  i  r 


PLOT    70 
for    applied    parameters    see   first    p 


page    of   this    appendix 
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PLOT    13 


ROLL 


q  kj  r 

1  •_ 


r 


uersus  n 


M 


x-src 


L 


.UL.        '-r  -±- 


UNITS    INCH., 


y-SCRLE:  2-0QE-01  UNITS    INCH* 

PLOT  71 
for  applied  parameters  see  first  page  of  this  appendix 
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W,'VMv. 


v'VVVWv^ 


K-3CPlLE=  i-OCF+Oi  UNITS  INCH. 
V-3CRL.F  -5.0CE.-0i  UNITS  INCH. 

RGRDTR    ,    TijRH   20   KH 
Pl3t    IS   RGll   RRTE 


RUD-15 

VERSUS   TIME 


PLOT    7  2 
for    applied    parameters    see    first    page    of    this   appendix 
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!  : 


1     nrr 


L. 


i  r? 


—  K 


T 


h     u 


■  r-        -i 


L.  u-  .        "  U 


•  N  L  ! 
llsiTT 


l  N  L 


l-i 


r-    \ 

J- 


PLOT    73 
for    applied    parameters    see   first    page    of    this    appendix 
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Ll, 


I      i       i      y         N<^.ft- 


'^r    v'S^^v'g 


1/  v 


Pl.^t    is  rqlL  rrte 


<  IFRSI 


r 


T  T  M 


V-SfC 


U 


/-1U 


rr     n~i 


Ui.       w  -L.       Ul  il    I  v        1.  1  1  s.;  1 


PLOT    74 
for    applied    parameters    see    first    page    of   this    appendix 
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Pi  FIT 


l_U 


T    ;- 
1   -^ 


ROLL 


r" 


:  irn c  i  i  o    t  t  k/i c 

'v      . —  1      !   'J  l»>  vJ  -L.I   ' 


PLOT    75 
for    applied    parameters    see   first    page    of    this   appendix 
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ft 

I 


I  (1, 


H-m 


!  iiv 


Pi  . 


HULL    r--       i 


IFR-- 


u  J 


T  I 


n  l 


5r^L; 


-  1 


U! 


r        t  i 


PLOT    76 
for    applied    parameters    see    first    page    of   this   appendix 
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■<■ 

-5CRLF-- 

1 

jE 

-Gi  UH 

r  75 

INCH, 

4 

-5CRLF.: 

^ 

n 
•  U 

-ip 

-01   J1-: 

i  !  j 

INCH* 

RGRO 

V 

'J 

r 

TUR 

H 

A  U    K  k 

P 

LOT 

I 

5 

RQL 

L 

r-  L  i  n ■  i     r 

RUQM-15 

VERSUS   TIME 


PLOT    77 
for    applied   parameters    see    first    page    of    this    appendix 
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\i  v'  w 


j:« 


K -SCALE ••■■!.  00F.+O1  JNITS  INCH, 
'-SCRLF.-G-OGE-Ql  JNITS  INCH 

RGR0K3    -    TURN   20   KN 
Pl'3T    13   RQlL   RRTE 


RUDM=15 


UFR 


c 


TTMF 


PLOT    78 
for    applied    parameters    see    first    page    of   this   appendix 
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V\Aa 


K-SCRLF.-  L.OCE+Qi  UNITS  INCH, 
V-3CRLE.-  2.0CE-01  UNITS  INCH, 

RGR1K3    r    TURN   2G   KH    ,    RUQM-15 
PlQT    IS        RHll   RNGLE  VERSUS   tIME 

PLOT  79 
for  applied  parameters  see  first  page  of  this  appendix 
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K~5rqLF.:\L03F+Qi  UNITS  INCH* 
t,-SnRLE----5-G0F.~0i  UNITS  INCH, 

RGR.1K3    -    TURN   20 
PLOT    15   ROlL   RRTE 


/:  k 


Hi   '  n  k/i        i    <~ 


!h 


JU 


D 


T 


PLOT    80 
for    applied    parameters    see   first    page    of    this   appendix 
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APPENDIX    B 


MODIFICATIONS    OF    THE    SIMULATION    PROGRAM 


-  During  the  course  of  studies  undertaken  with  the  XR-3 
Loads  and  Motion  Program  some  statements  have  been  added  or 
changed  in  order  to  improve  the  work,  with  the  simulation 
program. 

*  MAIN    program 

The  statements 

COMMON    /AIR/    PINF,RHOINF, GAM, IQUIT  (MAIN0050) 

IF     (MYTIME(DUM)  .LT. IQUIT)     GO    TO    12  (MAIN0661) 

have  been  modified  and  added  to  the  MAIN-program.  The 
purpose  of  the  second  statement  is  -co  start  the  writing  of 
the  output  values  desired  if  computer  time  has  reached  the 
demanded  time  (jobcard)  reduced  by  IQUIT  which  is  the 
expected  amount  of  time  (in  10  4  seconds)  required  to  print 
the  output.  The  user  may  input  the  desired  value  of  IQUIT  as 
the  fourth  parameter  (I 10-Format,  columns  36-45)  on  the  107 
data   card.       Default    value    for    IQUIT    is    600000    (60      seconds). 

*  Subroutine    INCON 

The  added  statements  are 

IQUIT=600000  (INCN0381) 

to    set   the   default    value   and 

READ     (5,2041)     IQUIT  (INCN1455) 
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to  read  the  set  value  for  IQUIT  thereby  overriding  the 
default  value. 

*  Subroutine  RHS 

In  Menzel's  version  [ Ref .  2]  of  the  XR-3  Loads  and  Motions 
Program  the  entered  craft's  velocity  00  (kn)  is  changed  in 
Subroutine  INCON  no  a  (ft/sec)  by 

U  =  00*1.6889  (INCN4580) 

and  in  Subroutine  RHS  transformed  back  to  units  in  knots 
(RHS  1850)  by 

VEL  =  0.5925*0  .  (RHS  1850) 

From    these    two   equations    follows 

VEL    =    0.592  5*00*1  .6889 

and  for  00=20  kn  there  results  VEL  =  20.0135  kn  .  But  since 
it  is  desired  that  VEL=00  the  statement  has  been  corrected 
to 

VEL    =    0/1 .6889     . 

Statements  MAIN0840,  RHS  1140  have  been  changed  accordingly. 


*  Subroutine  INTGRL 

To  this  subroutine  the  statement 

CALL  COLFIL  (INT  1021) 

has  been  added  in  order  to  provide  the  output  values  already 
calculated  to  the  user  if  the  minimum  time  interval  allowed 
in  the  integration  process  is  undergone.  Before  this 
addition  only  the  warning  message 

•DELTA  TIME  LESS  THAN  1.0E-6  -  -  JOB  STOPS' 

appeared  and  the  output  values  have  been  lost. 
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APPENDIX    C 


XR-3    LOADS    AND    MOTIONS    PROGRAM 
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